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19.  ABSTRACT  (Concluded) 


’piezoelectric  constant  e^?  and  the  elastic  constant  ,C12*  with  precise  measurements  of 
Che  oriencacion  of  crystals  to  be  used  in  the  experiments,  is  needed  before  drawing 
further  conclusions  about  the  causes  for  the  disagreement  between  theory  and 
experiment.  Specific  recommendations  for  future  theoretical  and  experimental 
efforts  are  suggested.  ^ _ 
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SUMMARY 


The  ability  to  make  dynamic  shear  stress  measurements  is  an 
important  need  in  both  laboratory  and  fielding  applications.  In  a 
preliminary  study,  the  feasibility  of  developing  a  piezoelectric  shear 
stress  gage  that  would  be  insensitive  to  compressive  stresses  was 
demonstrated.  The  objectives  of  the  present  work  were  to  perform 
detailed  analyses  to  understand  the  response  of  such  a  gage,  perform 
laboratory  calibration  experiments,  and  check  cross-axis  shear  effects, 
all  as  an  intermediate  step  for  developing  a  gage  package  for  field 
use.  Because  the  piezoelectric  constants  are  Important  in  determining 
the  optimal  gage  orientation,  a  related  objective  was  to  verify  the 
constants  reported  in  the  literature. 

To  meet  these  objectives,  theoretical  analyses  and  laboratory 
experiments  were  performed.  Impact  experiments,  designed  to  provide 
temporal  separation  between  the  compression  and  shear  response  of  the 
gage,  were  performed  on  163®-  and  165.5°-rotated  Y-cut  Lithium  Niobate 
(LiNbO^)  gages.  The  163®  Y-cut  orientation  was  chosen  because  it  is  the 
optimal  orientation  for  shear-stress-induced  response  based  on  the 
existing  set  of  piezoelectric  constants.  The  165.5°  Y-cut  orientation 
is  a  specific  direction  and  provides  a  uniaxial  strain  response  under 
compression-only  loading. 

The  experimental  results  showed  that  the  163°  Y-cut  orientation  is 
considerably  more  sensitive  to  shear  loading  than  to  compressive 
loading.  Also,  the  163®  Y-cut  orientation  appears  to  be  insensitive  to 
cross-axis  shear.  The  experimental  results  suggested  that  impact 
misalignment  (or  tilt)  and  crystalline  anisotropy  effects  need  to  be 
taken  into  account  in  analyzing  the  experimental  data.  As  expected,  the 
165.5®  Y-cut  orientation  showed  a  considerably  larger  response  to 
compression  than  did  the  163°  Y-cut  orientation.  The  large  current  rise 
that  occurs  following  the  initial  jump  in  current  suggests  the  need  for 
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modeling  the  effects  of  electromechanical  coupling.  Some  of  the  experi¬ 
mental  results  showed  scatter,  which  could  be  attributed  to  variations 
in  tilt  between  experiments,  minor  variations  in  the  orientations  of  the 
crystals,  or  both. 

In  analyzing  the  gage  response,  we  examined  only  the  initial  jump 
in  current  to  avoid  complexity  due  to  electromechanical  coupling. 
Quantitative  methods  were  developed  to  include  crystalline  aniostropy 
and  impact  tilt  contributions  because  these  were  determined  to  be 
significant.  Detailed  comparisons  between  the  experimental  results  and 
theoretical  predictions  were  attempted.  For  the  shear  loading,  the 
agreement  between  theory  and  experiment  is  good  (within  3  percent).  On 
the  other  hand,  for  compression  loading  theory  and  experiment  agree  only 
to  within  10  to  50  percent;  not  all  the  differences  have  been  resolved. 
Two  potential  sources  of  errors  that  need  to  be  examined  in  future 
studies  are  the  accuracy  limits  on  the  value  of  the  piezoelectric  con¬ 
stant  e22  and  the  uncertainty  in  the  orientation  of  the  rotated-cut 
crystals.  Also,  the  value  of  the  elastic  constant  needs  to  be 
established  accurately.  In  future  studies  an  accuracy  to  within  ±  0.1° 
is  required.  A  detailed  electrostatic  analysis  for  rotated-cut  crystals 
has  been  presented.  To  apply  this  analysis  to  experimental  data 
requires  solutions  for  the  two-dimensional  electrostatic  problem. 

Unless  this  analysis  is  carefully  applied  to  data  on  rotated-cut 
crystals,  interpretations  based  on  these  data  are  suspect.  Theoretical 
expressions  describing  the  effects  of  electromechanical  coupling  under 
combined  compression  and  shear  loading  are  derived;  however,  these 
expressions  were  not  used  in  analyzing  the  data.  This  latter  task 
requires  further  numerical  effort. 

The  present  work  has  provided  a  detailed  examination  of  the  shear 
gage  concept,  and  on  the  basis  of  the  work  described  in  this  report,  we 
believe  that  the  development  of  the  shear  stress  gage  is  a  realistic 
undertaking.  An  important  conclusion  is  that  improved  reconciliation 
between  theory  and  experiment  is  necessary  before  developing  a  field 
gage,  and  this  task  will  require  further  theoretical  analyses  and  more 
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precise  experimental  work.  Specific  steps  for  future  work  are  outlined 
in  this  report.  An  important  change  from  the  earlier  conclusions  of  the 
feasibility  study  is  the  redefinition  of  the  criteria  for  the  optimal 
gage  orientation.  The  sensitivity  of  the  gage  output  to  compression 
loading  varies  markedly  with  small  deviations  of  the  crystal  orientation 
from  the  optimum  orientation  for  measuring  shear  stress.  This  behavior, 
combined  with  practical  limitations  to  the  accuracy  with  which  the 
crystal  orientation  can  be  determined,  suggests  that  the  crystal 
orientation  to  be  used  in  a  field  gage  should  be  chosen  as  a  compromise 
between  minimizing  its  compression  sensitivity  and  accommodating  the 
practical  constraints  of  crystal  orientation  accuracy. 
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SECTION  1 


INTRODUCTION 

1 . 1  BACKGROUND 

The  measurement  of  dynamic  stresses  and  loads  is  an  important 
aspect  of  the  DNA  experimental  effort.  Many  different  types  of 
transducers  have  been  developed  and  are  in  use  in  both  laboratory  and 
field  programs.  Proceedings  of  a  recent  DNA  conference  describe  the 
status  of  instrumentation  currently  in  use.^  The  development  of  a  shear 
stress  transducer  for  dynamic  measurements  is  recognized  as  an  important 
need  by  workers  associated  with  dynamic  load  measurements. 

In  1978,  we  proposed  the  idea  of  developing  a  piezoelectric  shear 
stress  gage  that  would  be  insensitive  to  compressive  stresses.  A  small 
effort  was  undertaken  to  determine  the  feasibility  of  developing  such  a 
gage.2  Theoretical  conditions  for  developing  the  desired  gage  were 
formulated  and  quantified,  and  idealized  mathematical  relations  were 
derived  for  using  this  gage.  The  analytic  work  indicated  that  of  the 
various  orientations  in  alpha  quartz  and  lithium  niobate  (LiNbO^)  only 
the  163°-rotated  Y-cut  orientation*  appeared  to  be  suitable  for  use  as  a 
shear  gage.  Experiments  were  performed  to  confirm  that  the  gage  had  a 
high  sensitivity  to  shear  stresses  and  a  low  sensitivity  to  compressive 
stresses.  Although  this  work  was  mostly  qualitative,  it  demonstrated 
the  theoretical  and  experimental  feasiability  of  developing  a 
piezoelectric  shear  gage.  For  readers  not  familiar  with  the  feasibility 
study,2  a  synopsis  of  the  past  work  is  presented  in  Appendix  A.  A  brief 
background  pertinent  to  the  feasibility  study  is  also  included  in  this 
synopsis . 


*  1 
We  are  using  the  Institute  of  Radio  Engineers  (IRE)  nomenclature, 

which  is  recommended  for  piezoelectric  crystals.  Rotated  Y-cut 

orientation  denotes  rotation  about  the  crystallographic  X-axis. 


1.2  OBJECTIVES  AND  APPROACH 


The  objectives  of  the  present  effort*  were  to  perform  more  detailed 
theoretical  calculations  to  better  understand  the  gage  response,  perform 
laboratory  experiments  to  calibrate  both  the  compression  and  shear  res¬ 
ponse  of  the  gage,  check  cross-axis  shear  effects,  estimate  uncertain¬ 
ties  In  gage  measurements,  and  suggest  guidelines  for  developing  a  gage 
package  for  field  use.  The  present  work  was  envisioned  as  a  detailed 
examination  of  the  shear  gage  concept  and  an  intermediate  step  toward 
developing  a  field  gage  package.  Because  of  the  extensive  characteriza¬ 
tion  of  LiNbO^  both  in  ultrasonics^  ’  ■*  and  in  compressive  shock  wave 
experiments,^  no  serious  problems  were  anticipated.^  Although  the 
objectives  of  these  earlier  ultrasonic  and  shock  wave  measurements  were 
different  from  our  current  objectives,  we  felt  that  the  fundamental 
constants  of  LiNbO^  were  sufficiently  well  known  to  permit  a  straight¬ 
forward  calibration.  The  experimental  effort  was  to  provide  a  quantita¬ 
tive  check  for  our  initial  ideas.  As  discovered  during  our  work,  not 
all  of  fundamental  material  constants  (elastic  and  piezoelectric)  are 
known  with  sufficient  accuracy. 

For  our  laboratory  experiments,  we  chose  two  orientations:  163° 

Y-cut  and  165.5°  Y-cut.  The  163°  Y-cut  was  chosen  as  the  optimal 

direction  for  a  shear  stress  gage.  The  165.5°  Y-cut  orientation  was 

chosen  because  it  is  a  "specific  direction"^  and  theoretical  analysis  is 

easier  for  this  orientation.  The  165.5°  Y-cut  orientation  was  also 

intended  to  confirm  the  piezoelectric  constants  reported  In  the 
4-6 

literature.  We  felt  that  this  conformation  would  be  useful  because 


*Discussions  with  AFWL  personnel  were  helpful  in  determining  the  field 
requirements  for  a  shear  stress  gage. 

^As  described  in  Reference  6,  rotated-cut  gage  data  under  compression 
had  also  been  analyzed. 

"^A  direction  in  which  propagated  waves  are  purely  longitudinal  or 
purely  transverse  is  a  specific  direction.  Also  see  Appendix  B. 


the  optimal  orientation  for  the  shear  gage  is  sensitive  to  the  value  of 
these  constants. 


The  theoretical  work  was  intended  to  extend  the  analysis  presented 
in  our  feasibility  study  to  provide  sufficiently  accurate  quantitative 
interpretation  of  our  experimental  data.  A  secondary  objective  of  the 
theoretical  effort  was  to  Identify  additional  analyses  that  might  be 
important  for  the  subsequent  effort  of  developing  a  field  gage  package. 

The  approach  we  used  in  this  project  had  one  main  theme:  to  work 
out  as  many  of  the  quantitative  details  about  the  gage  response  as 
possible  to  minimize  ambiguities  and  difficulties  with  data 
interpretation  once  the  gage  is  used  in  field  experiments.  In  most 
applications  the  gage  will  be  subjected  to  complex  loading.  This  fact, 
coupled  with  the  fact  that  compressive  loads  are  considerably  higher 
than  shear  loads  in  most  applications,  requires  that  the  gage  response 
be  well  understood.  Hence,  it  is  necessary  that  in  well-defined 
laboratory  experiments  the  gage  response  be  modeled  quantitatively. 

The  remainder  of  the  report  consists  of  four  sections  and  several 
appendices.  Readers  desiring  an  overview  of  this  work  should  read 
Appendix  A,  the  first  few  pages  of  Section  2,  the  Summary  at  the  end  of 
Section  3,  Section  4.4,  and  Section  5.  The  experimental  work  is 
described  in  Section  3  and  Appendix  F.  Readers  interested  in  the 
theoretical  developments  and  analysis  should  read  Sections  2  and  4  and 
Appendices  B,  C,  D,  and  E. 


3 


SECTION  2 


THEORETICAL  DEVELOPMENTS 

In  this  section  we  present  theoretical  analyses  describing 
compression  and  shear  wave  propagation  in  piezoelectric  crystals  and 
develop  criteria  for  selecting  the  optimal  gage  design.  The  first 
theoretical  analysis  is  an  idealized  calculation  that  indicates  the 
requirements  for  an  ideal  shear  gage.  After  identifying  these 
requirements,  we  present  two,  more  general,  theoretical  analyses  that 
consider  conditions  that  are  closer  to  those  for  a  real  shear  gage. 
Although  these  two  analyses  have  not  been  fully  used  in  the  present  work 
because  of  their  complexity,  they  provide  considerable  insight  into  the 
gage  response.  Numerical  methods  to  Implement  these  general  analyses 
should  be  considered  in  future  studies. 

Since  Nielson's  initial  derivation  of  the  expressions  for  the 

short-circuit  current  for  a  quartz  gage,7  Graham  and  others  at  Sandla 

Laboratories  have  performed  an  extensive  series  of  studies  on  the 

response  of  piezoelectric  gages  subjected  to  dynamic  uniaxial  strain 
8—13 

loading.  Their  papers  demonstrate  the  extensive  developments  in 

both  theory  and  experiment  over  the  past  20  years.  In  contrast 

to  a-quartz  the  behavior  of  LiNbO^  is  considerably  more  complex  and 

requires  a  careful  understanding  of  the  electromechanical  coupling, 

which  is  discussed  in  Section  3  of  Reference  6.  The  papers  by  Chen  et 
14  15 

al.  and  Lawrence  and  Davison  are  important  contributions  to  this 
problem. 

Despite  the  many  developments  over  the  past  two  decades,  two 
features  are  common  to  all  the  studies  cited  above:  the  condition  of 
\niaxial  strain  and  the  assumption  of  coincidence  of  the  electric  fields 
with  the  gage  thickness  direction  (also  the  wave  propagation  direction). 
These  restrictions,  although  reasonable  for  the  earlier  work,  are  not 
valid  for  analyzing  the  shear  gage  response. 


*5 


All  the  theoretical  analyses  presented  here  assume  that  the  gage 
material  Is  elastic  and  has  zero  conductivity.  In  addition,  all  the 


dependent  variables  are  assumed  to  have  a  one-dimensional  variation. 
This  last  assumption  is  questionable  for  the  electrical  variables  in 
certain  situations,  as  discussed  later. 


2.1  ELASTIC-DIELECTRIC  ANALYSIS  FOR  COMBINED  COMPRESSION  AND  SHEAR 
LOADING 


Expressions  describing  the  response  of  an  idealized  gage  are 
obtained  by  extending  the  elastic-dielectric  analysis  for  uniaxial 
strain^  to  combined  compression  and  shear  loading.  We  make  the 
following  assumptions:  (1)  the  gage  thickness  direction  is  a  specific 
direction,  (2)  there  is  no  electromechanical  coupling,  (3)  finite  strain 
and  nonlinear  electrical  effects  can  be  neglected,  (4)  the  electrical 
permittivity  tensor  has  only  diagonal  elements*  and  their  variations 
with  strain  can  be  neglected,  and  (5)  the  electrical  fields  are  parallel 
to  the  gage  thickness  direction. 


With  these  assumptions,  we  consider  the  propagation  of  a 
compression  and  a  shear  wave  through  the  gage  as  shown  in  Figure  1. 
Region  I  is  unstressed,  Region  II  has  undergone  compressive  deformation, 
and  Region  III  has  undergone  compressive  and  shear  deformation.  The  two 
faces  of  the  gage  are  short-circuited  through  a  low  resistance  to 
measure  the  current  due  to  gage  deformation.  The  thicknesses  of  each  of 
these  regions  and  the  electrical  variables  for  the  three  regions  are 
indicated  in  Figure  1.  Because  there  is  no  free  charge  (i.e.,  the  gage 


has  zero  conductivity),  the  gradient  of  the  displacement  field  D  is 
_ 17 


(2-1) 


This  assumption  although  not  explicitly  stated,  has  been  made  in  past 
work. 


V  •'}  v->'- /W •V-.'*.  V'WVW  . 


where  X2  is  the  coordinate  through  the  gage  thickness;  the  gage  is 
assumed  to  be  infinite  along  the  other  two  directions.  Equation  (2-1) 
shows  that  D2  is  only  a  function  of  time  t,  and  at  any  instant  we  can 


write 


(1>D2  .  <2\  ■  <»„ 


(2-2) 


Writing  this  in  terms  of  polarization,  P,  and  electric  field,  E,  we  have 

d)£  (DE2  .  (2)e  (2)E2  +  C2)p2  .  (3)e  (3)E2  +  <3)p2  (2.3) 

where  is  the  permittivity  in  the  ic^  region.  The  short-circuit 

condition  is  given  by 


/ 1  E_dX.  -  ^E.A,  +  ^E-A„  +  ^E_A.  -  0 


'o  2  2 


(2-4) 


At  any  instant  the  thickness  A^,  A^,  A  can  be  expressed  in  terras  of  the 
total  thickness  A  and  the  wave  velocities: 


(1)E.(A  -  U  t)  +  <2>E„(U  -  U  )t  +  (3)E„U  t  -  0 

l  c  2  c  S  2  s 


(2-5) 


where  the  subscripts  "c"  and  "s"  refer  to  the  compression  and  the  shear 
wave,  respectively.  Because  we  are  considering  infinitesimal  strains, 
we  take  t  =  ^3^e  and  use  equation  (2-5)  with  equation  (2-3)  to 

eliminate  ^^E2  and  We  can  then  write 


<1>D2  ■  £  (1>E2  ’  T  '<2)p2(°S  -  V‘  +  <3>P2U,tl 


The  current  in  the  external  circuit  is  written  as 


1  =  A  dT  =  I  f(2)p^  +  u=(°)p9  -  (2)P->)1 


(2-6) 


2c  '  "s'  2 


(2-7) 
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Figure  1 .  Compression  and  shear  wave  fronts  in  a  piezoelectric  disk. 

[The  wave  velocities  Uc  and  U$  and  the  electrical  variables 
in  the  three  regions  are  shown.  The  gage  thickness  direction 
is  X2  (or  X2)  as  explained  in  the  text.  The  polarization 
in  region  I  is  zero.] 
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where  A  is  the  area  of  the  electrode.  In  equation  (2-7),  v  'V*  is  the 


polarization  due  to  the  compressive  wave  only,  but 


(3)f 


2 

is  the 


polarization  due  to  both  waves.  Rewriting  equation  (2-7),  we  have 


i  ■  7  [P  U  +  P  U  ] 
l  1  c  c  s  SJ 


(2-8) 


The  polarization  is  related  to  the  mechanical  strains  (in  the  gage) 
through  the  relation^-® 


6ijk  ejk 


(2-9) 


If  X2  is  the  wave  propagation  direction  and  the  direction  of  shear 
loading  is  X3,  we  can  rewrite  equation  (2-9),  in  the  matrix  notation,  in 
Region  III  as 


(3), 
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c  s 


e22  e2  +  e24 


(2-10) 


In  equation  (2-10)  and  throughout  the  report,  we  follow  the  commonly 
used  convention  for  equivalence  of  the  subscripts  in  matrix  and  tensor 
notation: 
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Substituting  equation  (2-10)  into  equation  (2-8),  we  can  write 


c  22'2 


s  24  4- 


(2-11) 
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Equation  (2-11)  can  be  used  to  discuss  the  criteria  for  selecting  a  gage 
subjected  to  combined  compression  and  shear  loading.  For  a  gage 
sensitive  only  to  normal  strains,  should  be  zero.  In  contrast,  an 
ideal  shear  gage  should  have  ^22  *  0*  To  satisfy  these  conditions,  we 
need  to  consider  rotated  cuts  for  the  crystals. 

We  choose  the  following  notation  to  avoid  confusion  with  the  axes 

» 

used  in  the  crystallographic  system:  X2  defines  the  direction  of  wave 
»  ♦ 
propagation,  X3  is  the  direction  of  the  imposed  shear  motion,  and  X^  is 

I 

chosen  such  that  X^-axes  form  a  right-handed  coordinate  system.  Hence, 

f  f 

an  important  requirement  for  an  ideal  shear  gage  is  that  e22  *  0  and  e?^ 

I  I 

be  large.  For  the  most  general  case,  we  have  shear  along  the  X^-  and  X3 
axes.  For  an  ideal  gage  there  should  be  no  cross-axis  shear  effect  and 

t 

e26  *  0*  For  the  situation  shown  in  Figure  1,  the  current  output  can 
then  be  uniquely  related  to  the  shear  strain  in  the  gage.  The  use  of 
this  idealized  gage  would  be  analogous  to  the  piezoelectric  compression 
gage. 


Having  established  the  main  criterion  for  the  shear  gage,  we  now 

examine  the  various  assumptions  made  in  deriving  equation  (2-11).  Of 

the  five  assumptions  listed  earlier,  the  most  important  one  is  the  need 

19 

for  the  wave  propagation  direction  to  be  a  specific  direction  with 

regard  to  both  compression  and  shear  waves  (see  Appendix  5).  If  the 

direction  of  wave  propagation  is  not  a  specific  direction,  the  various 

strains  (or  stresses)  in  the  gage  will  be  coupled,  all  three  strains 

(€',  e! ,  el)  will  be  produced,  and  the  polarization  can  have  several 
4  4  o 

contributions. 


Thus,  the  two  most  Important  criteria  for  the  shear  gage  are  as 

follows:  (1)  the  electrical  polarization  is  caused  only  by  the  shear 

strain  (or  stress)  of  interest  and  (2)  the  gage  thickness  direction  is 

2 

specific  direction.  In  the  earlier  feasibility  study,  these  criteria 
were  used  to  examine  both  a-quartz  and  lithium  niobate  crystals.  Only 
the  163°  Y-cut  orientation  of  LiXbO^  approached  optimal  conditions  (see 
Appendix  D).  For  this  orientation,  using  the  constants  cited  in  the 
literature,  we  find  that  the  ratio  ^2k^Q’22  *s  ’  Although  this 

9 


V  ■< 


«\‘»v 


direction  is  not  a  specific  direction,  it  is  close  to  the  specific 
direction  (165.5°  Y-cut,  see  Appendix  B).  The  particle  motion,  for  wave 
propagation  along  the  163°  Y-cut  in  LiNbO^,  is  along  the  161.7°  Y-cut 
direction.  Hence,  the  mechanical  coupling  of  the  strains  due  to  crystal 
anisotropy  was  expected  to  be  small. 

Note  that  the  values  cited  above  were  based  on  the  values  for 
elastic  constants  (clj’j)  and  piezoelectric  constants  for  LiNbO^ 

cited  in  the  literature.  The  extreme  sensitivity  of  the  ratio  ®24^e22 
to  the  crystal  orientation*  requires  that  both  the  ejj^  values  and  the 
crystal  orientation  be  known  very  accurately.  During  the  work  described 
in  this  report,  we  discovered  that  the  values  cited  in  the  literature 
were  incorrect.  We  have  reanalyzed  past  shock  wave  data  for  LiNbO^  to 
determine  the  "correct"  values  (see  Appendix  C).^ 

In  the  present  work  we  also  found  that  small  deviations  from  the 
specific  direction  and  the  effects  of  impact  tilt  need  to  be  included 
for  analysis  of  experimental  data.  Analytic  methods  to  incorporate 
these  effects  are  described  in  Appendices  B  and  E. 

In  the  next  two  subsections,  we  present  more  general  theoretical 
developments  that  are  closer  to  the  conditions  of  the  experiments. 

2.2  ELASTIC-D I ELECTRIC  ANALYSIS  FOR  A  ROTATED-CUT  GAGE 

In  this  subsection,  we  restrict  ourselves  to  only  two 
assumptions:  (l)  electromechanical  coupling  is  not  considered  and  (2) 

finite  strain  and  nonlinear  electrical  effects  are  neglected.  The  first 
assumption  is  relaxed  in  the  next  subsection.  Eliminating  the  second 
assumption  does  not  pose  any  conceptual  problems,  but  makes  the  analysis 


it  »  »  i  » 

The  ratio  624^e22  ^or  *62°  Y-cut  is  -28#8  and  the  ratio  Goh'^22  ^or 
164°  Y-cut  is  25.9.  Note,  is  nearly  constant,  but  the  e?7  value 

changes  dramatically  with  orientation. 

^This  question  has  not  been  resolved  satisfactorily. 


more  complicated.  For  the  small  strains  considered  in  our  work,  the 
second  assumption  is  reasonable. 

Because  we  are  not  considering  wave  propagation  in  a  specific 
direction,  each  of  the  Regions  II  and  III  in  Figure  1  will  have  both 
compressive  and  shear  deformations.  In  general,  there  will  be  three 
regions;  however,  to  simplify  the  algebra,  we  initially  consider  only 
two  regions:  an  unstressed  region,  Region  I,  and  a  region  undergoing 
both  compressive  and  shear  deformation,  Region  II.  Thus,  Region  II  in 
this  analysis  is  similar  to  Region  III  of  Figure  1.  Extending  this 
analysis  to  three  regions  is  straightforward  and  the  final  expressions 
are  presented.  The  two-region  situation  arises  when  a  quasi¬ 
longitudinal  wave  propagates  in  the  crystal;  this  situation  is 
encountered  in  many  of  our  experiments  described  in  Section  3.  Below  we 
derive  the  equations  for  a  rotated  Y-cut  crystal  of  LiNbOj  with  the 

f 

X2~axis  defined  to  be  the  wave  propagation  direction  in  Figure  1. 

Proceeding  as  before,  the  absenc  of  free  charge  in  the  crystal 
provides 


(2-12) 


Thus,  D2  is  only  a  function  of  time.  Writing  equation  (2-12)  in  terms 
of  $  and  £  does  not  produce  equation  (2-3)  as  before  because  for  rotated 
Y-cut  LiNbO^  crystals  the  permittivity  tensor  has  the  form 
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where  0  Is  the  angle  of  rotation  about  the  X^-axis.  It  Is  the  non- 
diagonal  form  of  the  permittivity  tensor  that  Introduces  most  of  the 
complexities  In  the  analysis  of  these  rotated-cut  crystals.  Using  the 
form  of  e  In  equation  (2-13),  we  can  rewrite  equation  (2-12)  as 


(Dp,  .  (1)  .  (DF,  ,  (ID,  (Dp,  .  (2)  ,  (2)  (2) 

P2  +  e22  E2  +  £23  E3  P2  +  e22  E2  (2-14) 


+  (2)  (2) 

e23  *"3 
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In  the  unstressed  region,  the  polarization  v  is  zero.  In  accordance 

with  assumption  (2)  above,  we  neglect  changes  in  the  permittivity  due  to 
small  strains.  The  continuity  of  the  tangential  component  of  the  E 
field  across  the  wave  front  gives 


(2-15) 


Equation  (2-14)  can  now  be  solved  to  give 


(2-16) 


The  short-circuit  condition  across  the  electrodes  on  the  gage  faces  is 
written  as 


e;  <“2 


,  a  (D-,  5 ,  ,  (2 ) „ ,  ,  n 

“*  t.  X.,  T  X..  =  u 


'2  1 


'2  2 


Solving  equations  (2-16)  and  (2-17)  gives 

(2)p,  ,, 

(Dr,  _  2  *2 

e2  vTT  ■ 


'22 


(2)p,  ., 
(Dp,  „  D  _1 

2  ' E  2  2 


(2-17) 


(2-18a) 


(2-18b) 
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In  the  small  strain  approximation,  we  are  neglecting  particle  velocity 

?  I 

uj  in  comparison  with  the  wave  velocity  Uc,  and  we  can  write 


V  =  l'  -  U't 
1  c 


and 


a:  =  u't 
2  c 


(2-19) 


Here  the  subscript  "c"  refers  to  the  quasi-longitudinal  wave.  From 
equation  (2-17)  we  see  that  the  electric  fields  depend  on  the  external 
circuit  condition;  equation  (2-19)  shows  that  they  are  strongly  time- 
dependent.  Combining  equations  (2-1P)  and  (2-19),  we  can  write 


^P' U't 

ay _ X£-,£.  ay 

2  l'  e23  j 


(2-20) 


Note,  in  contrast  to  equation  (2-6),  the  displacement  field  in  the  gage 

I 

thickness  direction  depends  on  the  electric  field  along  the  X^- 
direction.  The  current  in  the  external  circuit  is  the  displacement 
current,  which  is  given  by 


i  =  A 


d°2 


dt 


—  r a ) p i  i  +£'e*  a y ' 

r  1  r2  c  23 


(2-21) 


where  the  dot  denotes  a  time  derivative.  We  write  the  expression  for 
current  as 


(2)p;au’ 


(2-22) 


( I)  • 

where  the  current  i_  =  A  e'  E'.  For  a  rotated  Y-cut  LiNbO^  crystal 

t  23  3  J 


the  onlv  nonzero  strain  components  are  e',  e' ,  and  e' .  Hence,  the 

.24  6 


expression  for  polarization  along  the  X2~direction  is  given  by 


(2) 

V  'V<  -  0>  C  <  +  ='  P<  +  a’ 

2  22  2  24  4  26  6 


(2-23) 


It  can  be  shown  that,  for  a  rotated  Y-cut  LiNbO^  crystal,  the  constant 

I 

e26  =  Hence,  the  last  term  on  the  right  side  of  equation  (2-23)  is 

zero.  Substituting  equation  (2-23)  in  (2-22)  produces 


From  equation  (2-24)  we  see  that,  in  the  approximation  of  no  electro¬ 
mechanical  coupling,  except  for  ig,  the  current  produced  by  wave 
propagation  in  a  nonspecific  direction  is  similar  in  form  to  the 
previous  derivation  in  Section  2.2.  The  simpler  relationship  holds  for 

Y-cut  LiNbO-,  crystals;  the  term  i_  is  zero  because  e*  only  has  diagonal 
j  t,  ij 


terms . 


( 2 )  • 

To  determine  the  effect  of  E^  on  the  current  in  equation  (2-21) 

is  an  involved  problem.*  It  will  depend  on  the  electrical  boundary 

» 

condition  along  the  X^-direction.  For  a  truly  infinite  disk,  this  term 
will  not  contribute.  However,  for  a  finite  disk,  the  shunted  and 
shorted  gages  may  show  markedly  different  behavior.  Further  analysis  of 
this  problem  requires  a  two-dimensional  electrostatic  analysis  and  is 
not  attempted  here.  In  future  work  this  issue  needs  to  be  carefully 
examined. 


For  Y-cut  and  rotated  Y-cut  LINbO-,  crystals  there  is  an  additional 
(  2  )  • 

factor  influencing  El  and  i-.  For  these  orientations  there  is  also  a 

^  i 

nonzero  polarization  in  the  X^-direction,  even  for  compression-only 

loading.  This  polarization  may  alter  the  electrical  boundary 

condition.  Thus  Its  magnitude  is  also  needed  to  solve  the  electrostatic 

C  2 )  • 

problem  to  determine  . 

If  we  consider  three  regions^  as  shown  in  Figure  1,  we  can  write 


ll 


l'  -  U't 
c 


l ' 
2 


(U'  -  U')t 
c  s 


U't 

s 


(2-25) 


★ 

Discussions  with  Mr.  R.  A.  Graham  of  Sandia  National  Laboratory  on  this 
problem  have  been  helpful. 

+ 


The  most  general  case  would  have  four  regions  because  there  will  be 
two  quasi-shear  waves.  For  the  present  problem,  one  quasi-shear  wave 
is  sufficient. 


i 


► 


K 


Here,  Uc  and  Ug  are  the  wave  velocities  of  the  quasi-longitudinal  and 
quasi-shear  wave,  respectively.  Proceeding  as  before,  the  displacement 
field  can  be  written  as 


(1)D^  «  jr  [<2)P£  (IT  -  U')t  +  (3)P£  U’t]  +  e^3  (1)E^  (2-26) 

Comparing  equations  (2-6)  and  (2-26)  we  see  that  the  two  are  similar 
» 

except  for  the  E^  term.  The  current  in  the  external  circuit  can  be 
written  as 


dD£ 

dt 


—  { 

V  1 


(2) 


.(3) 


(2) 


P2Uc  +  *  P2  "  +  re23  EP 


(2-27) 


—  rp '  u’  +  p*  n*  i  +  i 

V  1  2c  c  *2s  UsJ  E 


The  values  of  i^  in  equations  (2-22)  and  (2-27)  will  be  different 
because  the  polarization  fields  in  the  X-j -direction  will  be  different  in 
the  two  cases. 

2.3  LINEAR  COUPLED  RESPONSE  OF  A  ROTATED-CUT  GAGE  TO  COMBINED 

COMPRESSION  AND  SHEAR  LOADING 

We  now  present  the  theoretical  formalism  for  including  electro¬ 
mechanical  coupling  in  the  gage  output  analysis.  The  theoretical 
developments  for  compression-only  loading  and  without  the  complexity 
introduced  by  rotated  cuts  can  be  found  in  the  paper  by  Chen  et  al.*^ 
However,  these  authors  include  finite  strain  and  more  general  external 
electrical  circuit  conditions  in  their  work. 

The  experimental  situation  analyzed  here  is  comparable  to  that 
shown  in  Figure  l.  The  approach  to  solving  the  coupled  problem  is 
similar  to  that  for  other  wave  propagation  problems.  The  constitutive 
equation  in  combination  with  the  governing  equations  yields  a  wave 
equation,  which  is  then  solved  to  provide  the  wave  profiles  at  different 
positions  in  the  sample.  Except  for  a  linear  elastic  solid,  a  numerical 


method  is  needed  to  obtain  the  wave  profiles.  In  the  coupled  problem, 

the  constitutive  equation  has  both  mechanical  and  electrical  terms. 

Hence,  the  external  circuit  conditions  can  influence  the  wave  profiles. 

A  comprehensive  discussion  of  this  topic  is  presented  in  an  article  by 
20 

Thurston. 

For  the  one-dimensional  problem  of  interest,  we  can  write  the 
20 

governing  equations  as 


(2-28) 


(2-29) 


(2-30) 


(2-31) 


f  t 

Here  0q  is  the  initial  density,  u^  is  the  particle  velocity  in  the  X^- 
direction,  and  £  is  the  specific  internal  energy.  We  are  using  the 
small  strain  approximation  and  considering  stresses  and  strains  to  be 
positive  in  tension. 


By  choosing  the  strain  (e^)  and  the  displacement  field  (Dm) 
Independent  variables,  we  can  write  the  constitutive  equation  for 
linear  piezoelectric  solid  in  the  matrix  notation  [see  the  matrix 
tensor  notation  following  equation  (2-10)]* 


as 

a 

and 


*The  constitutive  equations  are  in  the  Xj  system  unlike  the  governing 
equations  (2-28)  through  (2-31). 


16 


a 

a 


-  h  D 
mat  m 


(2-32) 


E 


i 


€  +  Bf. 

a  ij 


D  . 
J 


(2-33) 


where 


and  are  the  elastic  stiffness 
assume  isentropic  deformation  and 
can  be  written  as 


constants  in  a  constant  $  field.  We 
recognize  that  the  Maxwell  relations 


(2-35) 


Further  details  about  the  constitutive  equations  can  be  found  in  Chapter 
4  of  Mason's  book^  or  Section  14  of  Thurston's  review  article.^®  The 
constitutive  equations  clearly  show  the  dependence  of  the  elastic  moduli 
and,  hence,  the  wave  velocities  on  the  electrical  variables.  In 
Appendix  B,  the  relation  between  C®j  and  c![j  is  derived. 


In  the  absence  of  electromechanical  coupling,  equations  (2-28) 
through  (2-31)  in  conjunction  with  the  constitutive  equations  can  be 
solved  for  a  particular  initial  condition  using  a  generalized  one¬ 
dimensional  wave  propagation  program. ^  Such  calculations  are  described 
in  Appendix  B.  For  the  electromechanical  problem  we  need  to  specify  the 
external  electrical  conditions.  A  simple  resistive  circuit  is 
considered. 


The  current  flowing  in  the  external  circuit  is  given  by 


i(t) 


(2-36) 


The  voltage,  across  a  resistor  R,  in  the  external  circuit  is 


17 


V(t)  -  ARD,  -  -  /  e:  dX 


2  2 


(2-37  ) 


For  a  rotated-cut  crystal,  we  have 


pi  +  eIj  Ej 


(2-38) 


We  can  eliminate  E2  between  equations  (2-37)  and  (2-38)  and  write  a 
differential  equation  for  D2(t).  We  consider  the  specific  situation  of 
a  rotated  Y-cut  LiNbO-j  crystal  and  use  the  permittivity  tensor  from 
equation  (2-13).  Combining  equations  (2-37)  and  (2-38),  the  voltage  in 


the  external  circuit  is 


V  Di  -  Pi  -  el,  El 

-  /  - — dx; 

0  E22 


(2-39) 


Because  of  equation  (2-12),  D2  is  only  a  function  of  time.  The 
tangential  component  of  the  E  field  is  continuous  through  the  thickness 
of  the  disk.  As  before,  we  assume  the  permittivity  to  be  unchanged  by 
the  small  strain.  Hence, 


(°2  *  e23  E3>*’ 


P2  dX2 


(2-40) 


P,  in  the  disk  depends  on  all  three  strain  components  (e',  e ' ,  and  e') 
as  indicated  in  equation  (2-23). 

The  general  solution  of  the  coupled  problem  is  obtained  by  solving 
equations  (2-28)  through  (2-30),  (2-32),  and  (2-40)  simultaneously. 
Lawrence  and  Davidson^  presented  a  numerical  method  for  solving  the 
coupled  problem  for  the  simpler  uniaxial  strain  loading.  This  solution 
will  be  useful  as  a  guide  in  developing  the  solution  for  the  more 
complex  situation  considered  here. 


The  form  of  this  equation  is  identical  to  Ea.  (3.2)  in  the  paper  by 
Chen  et  al.  in  the  limit  of  small  strains.  **  A  similar  equation  was 
derived  20  years  ago  by  Jones  et  al. 


Cl? 

w'KlS-.)' 


In  the  short-circuit  approximation  (R  *  0),  we  obtain  the  results 
derived  in  the  last  subsection. 


°2  /0  P2  ^2  +  e23  E3 


(2-41) 


Although  this  expression  is  the  same  as  that  derived  for  the  elastic- 
dielectric  medium  in  the  uncoupled  approximation,  the  similarity  is 
somewhat  misleading  because  the  values  of  the  D  field  at  any  instant 
determine  the  mechanical  variables  through  equation  (2-32).  Hence,  the 
electrical  and  mechanical  variables  need  to  be  calculated  simultaneous¬ 
ly.  Figure  4  in  Reference  15  shows  clearly  the  effect  of  electro¬ 
mechanical  coupling  on  the  current  output  from  a  Z-cut  LiNbC^  gage.  The 
electrical  field  variations  can  produce  stresses  in  regions  where  the 
mechanical  wave  has  not  reached.  To  include  electromechanical  coupling 
in  the  analysis  of  wave  propagation  (presented  in  Appendix  B)  in 
rotated-cut  crystals  is  a  time-consuming  but  straightforward  procedure. 
However,  before  such  a  task  is  attempted,  the  importance  of  the  two- 
dimensional  electrical  field  effects  indicated  at  the  end  of  the  last 
subsection  needs  to  be  determined. 

The  electromechanical  coupling  effects  can  be  neglected  in  two 
situations:  in  the  analysis  of  the  initial  current  jump  as  discussed  by 

Graham*’  and  in  the  open-circuit  mode  where  no  current  is  allowed  to  flow 
In  the  external  circuit.  This  latter  situation  arises  in  the  field 

usage  of  piezoelectric  gages.  In  the  open-circuit  mode,  i  *  0  in 

»  * 

equation  (2-36).  Hence,  D2  is  independent  of  time.  Thus,  the 
mechanical  variables  can  be  calculated  in  the  usual  manner  by  using 
as  the  elastic  moduli.  For  the  open-circuit  case,  the  mechanical 
boundary  conditions  and  their  effects  on  the  electrical  fields  may  be 
important . 


Because  of  zero  conductivity,  D2  Is  independent  of  Xj ■ 


2.4  IMPLICATIONS  FOR  EXPERIMENTS 


It  is  desirable  to  find  a  gage  whose  output  can  be  modeled  using 
the  idealized  elastic  dielectric  analysis  presented  in  Section  2.2. 

Such  a  gage  would  have  to  meet  the  following  three  requirements:  (1) 
its  electrical  polarization  must  be  due  solely  to  the  shear  stress  of 
interest  (no  cross-axis  shear  effect),  (2)  this  polarization  must  be 
one-dimensional  and  along  the  gage  thickness  direction,  and  (3)  the  gage 
thickness  direction  must  be  a  specific  direction. 

2 

As  indicated  in  our  previous  study,  none  of  the  principal  cuts  in 
a-quartz  or  LiNbf^  satisfy  all  of  these  requirements.  Even  when 
rotated-cuts  are  used*  only  the  163“  Y-cut  LiNbO^  crystal  approaches  the 
desired  requirements.  As  indicated  earlier,  the  optimal  orientation 
depends  strongly  on  the  values  of  the  piezoelectric  constants  e^j. 

Also,  for  the  gage  to  be  useful  the  crystal  orientation  has  to  be 
matched  accurately  to  its  optimal  value.  Practical  usage  always 
involves  some  error  or  variation.  Hence,  laboratory  experiments  are 
required  to  assess  the  gage  response,  to  verify  the  theoretically 
predicted  behavior,  and  to  determine  the  tolerance  limits  for  the  gage 
response.  Because  the  output  of  the  rotated-cut  crystals  is  not  easily 
analyzed,  the  experiments  are  designed  using  the  Idealized  analysis. 
Subsequently,  the  gage  output  is  carefully  analyzed  to  assess  the 
deviation  from  the  idealized  behavior. 

Although  the  current  or  short-circuit  mode  is  more  complex,  it  is 
ideally  suited  for  laboratory  calibration  studies-  The  initial  jump  can 
be  analyzed  without  the  complexity  of  electromechanical  coupling.  The 
results  can  then  be  used  in  developing  a  gage  package  for  field  applica¬ 
tions. 

Our  feasibility  study  indicated  that  the  specific  direction  closest 

2 

to  the  163°  Y-cut  orientation  was  the  165.5°  Y-cut  orientation.  The 


*The  accepted  nomenclature  for  piezoelectric  crystals  is  presented  in 
Reference  3. 


closeness  of  these  two  orientations  suggested  that  they  both  should  be 
examined  in  the  laboratory  experiments.  The  one-dimensional  compression 
and  shear  experiments  we  performed  on  gages  with  these  orientations  are 
described  in  the  next  section.  These  experiments  allowed  us  to  quanti¬ 
tatively  measure  the  shear-to-compression  sensitivity  ratio  (163°  Y-cut 
orientation),  check  the  piezoelectric  constants  by  performing  experiments 
along  a  specific  direction  (165.5°  Y-cut  orientation),  and  determine  the 
effect  of  small  deviations  from  a  specific  direction.  In  addition, 
detailed  analysis  of  the  gage  response  can  be  used  to  ascertain  the  need 
to  model  the  more  complex  phenomena  indicated  in  this  section. 


.V 


SECTION  3 


EXPERIMENTAL  METHOD  AND  RESULTS* 

3.1  EXPERIMENTAL  DESIGN  AND  METHOD 

All  the  LiNbO-j  gages  used  In  our  work  were  fabricated  by  Specialty 

Engineering  Associates  according  to  our  specification.  The  LiNbO^ 

crystals  were  cut  to  the  desired  orientation  to  within  ±  0.5°.^  The 

gages  were  parallel  to  within  ±  2.5  pm  and  were  planar  to  within 

±  1  pm.  These  values  are  typical  of  piezoelectric  gages  used  In  shock 

wave  experiments.  The  wraparound  (or  shorted  gage)  configuration  was 

used  In  all  the  experiments  described  here.  A  few  experiments  were 

performed  with  the  shunted-gage  configuration,  but  It  was  difficult  to 

maintain  a  ground  connection  with  the  gage  front  surface.  Because  the 

23 

shunted  gage  configuration  Is  a  superior  configuration,  the 
possibility  of  using  shunted  gages  in  future  experiments  should  be 
examined. 

Guard  rings,  nominally  0.08  mm  wide,  were  cut  on  the  back  side,  and 
a  50-Q  resistor  was  connected  across  the  inner  electrode  and  the  ground 
(or  outer)  electrode.  The  gages  were  enclosed  in  an  aluminum  housing 
with  epoxy  potting  at  the  rear  of  the  gage.  The  nominal  values  for  gage 
diameter,  gage  thickness,  and  the  inner  electrode  diameter  were  26  mm, 

5  mm,  and  8  mm,  respectively.  The  exact  values  for  the  different  gages 
along  with  the  experiment  number  are  presented  in  Table  1.  Some  of  the 
columns  in  Table  1  are  explained  in  the  following  paragraphs. 


The  work  described  in  this  section  was  performed  in  collaboration 
with  W.  J.  Murri. 

^The  supplier  indicated  that  the  orientations  were  good  to  within 
±  0.2°.  We  now  believe  that  future  work  will  require  that  the 
orientations  be  known  with  better  accuracy. 


? 

i 

i 


* 

\ 

I 

F 

r 

r 

l 

r 

l 

I 

r 

t 

i 


*  pnw  win  uLr^wya-i  wmi -  -*n  i 


*•’*-0  w-'-WTri 


! 

V 

04 

CM  CM  CM  OJ 

04  04  04  04  04  04 

04 

04  04 

«l 

> 

V  o 

o 

o  o 

^  —4 

o  o 

o  o 

«— 4  *-4 

o  o 

o  o 

P-*  4-4 

© 

r-4 

o  © 

H 

VI 

fi 

>  ns 
<004 

X 

X  X 

X  X 

X  X 

X  X 

X  X 

X 

X  X 

Jrf 

y 

3  9 

av  so 

os  OS 
n  co 

<n  x 

CO  CO 

m  os 

n  n 

^  04 

y  -y 

CO 

y 

in  os 
rs  rs 

X 

y  -y 

-y  <r 

«y  -y 

•y  «y 

y  y 

f- 

<0 

m 

J,  Jt 

J,  ^ 

»— 4  9-4 

>—4  —4 

0-4 

k. 

«< 

x 

V 

y 

in 

in 

X  in 

s*« 

V 

rs  oo 

O 

-y  O 

O  o> 

O  O' 

04 

—  © 

44 

y 

V 

44 

09  a 

-4  w 

4A 

CSl 

<T  X 

CM  04 

\0  n 

04  OJ 

x  x 

04  04 

x  x 

04  04 

in  n 

Ol  04 

X 

04 

X  X 

Ol  Ol 

<4* 

00 

? 

V 

td 

00 

1 

4J 

c 

V  <o 

u  ; 

in 

x  in 

a  u 

V 

4-1  1 

o-  O 

o*  04 

04  -y 

O'  x 

oJ  y 

X 

os  rs. 

U  V 

(X 

«0 

o 

)  ®  nv 

•3  3 

O 

o  rs 
n  n 

o  04 
n  n 

-2  o 

n  n 

Os  -- 
-y  in 

04  04 

in  in 

04 

in 

CO  CM 
in  in 

V  X 

a.  n 

V 

l 

V  ^  o 

00 

u  *-» 

1 

44  <0 

g  41  b 
to  X  <0 

at 

9 

a  «  ^ 

s/i 

•** 

* 

X 

4-1  4-S 

■*3  g 

09 

O 

-<  O 

0O  x 
o  o 

os  -» 
o.  X 

o  o 

O' 

-  o 

^  s o 
oo  X 

o  o 

X  OS 

X  O' 

°  °. 

O' 

© 

y 

O'  O' 

o  © 

e  s 

s*  14*  O 

0  -o 

V  X 

30 

T3 

U 

1  -*  9 
3  w- 

© 

o’  o 

o  o 

o  o 

o  o 

o  © 

© 

©  © 

x  c  e 

44  O  V 

S3 

« 

1 

a 

50 

2 

1 

C  44  U 

o 

W  O  V 

v  a 

o 

, 

Qfi  U  s* 

2 

V  u 

C  s-l 

z 

t: 
u  c 

V  k 

5  V 

>  44  ^ 

4  v  a 

8 

x  ai 
<3>  ai 

X  -N 

O'  O 

8  S 

O'  X 

Os  O' 

Q  Os 

O  O' 

O' 

O' 

88 

^4  x  c 
a  v  o 

□ 

C  4 

c  C, 

J  B  1 

l  as  W 

-T 

c-» 

o.  X 

00  co 

os 

OO  Os 

Os 

X  X 

o  c  u 

x  o  y  • 

>-<  V  -4 

s-4  OJ  V 

-*  a 

IJ  4-4  U  CL 

Vi 

CxJ 

M  U  -  o 

<tj  v  x  u 

00  b  <0 

o> 

s-  0 

© 

c 

CO 

| 

V  X  <0  o- 

V 

«  c  ^ 
oO-x  g 

S 

es  <m 
o  o 

04  O' 

O  <ys 

X  O' 

O'  O' 

04  04 

O  O 

O'  Os 

O'  O' 

O' 

O' 

X  X 

OS  O' 

X  -- 

4J  W  C 

to  o  y 

, 

<0  0  9 
O 

m  n 

n  n 

n  ~y 

■y  -y 

in  n 

y  y 

y 

y  y 

c  y  n  co  • 

O  -4  --  O  <N 

«■» 

X 

| 

w  u 

H 

1 

«  W  0  V  c 

ss  V  X  X  0 

X  >  4J  — 

^5 

H 

u 

<0  V  44 

V  * 

V  w  nsl  O 
IM  II  O 

«n  n 
o  o 

n  -y 
O  O 

~y  -t 

o  o 

n  n 

O  O 

y  y 
©  © 

y 

© 

y  y 

O  © 

f  V  X  4*  o 

-  X  4J  <0  V 

N  4-*  X 

t 

SB  HI 

>  *  w:  S 

X  x 

CM  04 

45  X 
04  04 

X  X 
04  04 

X  X 
04  04 

x  x 

04  04 

X 

04 

X  X 

04  04 

00  c 

VO0CO**- 
X  c  o  -o  o 

u 

1 

*J  O  4-1 

0*  <0  <0  <^N 

14*  <0  3  X 

O  «  —  I 

<0 

■ 

v>  <0  e  04 

Ci 

c  <0  3  ^ 

O! 

0  3  « 

u  —  1 

«  Ji  e 

X  <0  -  O 
<  *■*! 

0  o 

O  O 

9  0 

o  o 

9  9 

o  o 

9  0 

o  o 

Os 

9  9 

©  © 

O' 

9 

o 

O' 

9  9 

o  © 

O' 

—  DUG 

4-4  'J3  -H  0 

4J  X  <0  X  *J 

1  c 

C  <0  1  *-•  <0 

-  v . 

V  1  -  3 

09  —  > 

—  -  CM  4-i  CT 

U  1 

U  CsJ  C  V 

C 

51 

0  V  3 

V  X  0  c 

v  X  Ol  U  — 

cut 

cut 

cut 

cut 

Cut 

W  4-1 

3  3 

<J  <J 

u  u 
3  3 

<J  <J 

4J  4-4 

3  3 

(J  o 

4-4 

3 

O 

4-4  4-4 

3  3 

U  V 

X  44  O 

4-4  w  -0 

si  "0  V 

o  !0  x  0  C 

V  v  1  > 

>-  > 

>•  >- 

>-  >« 

>4  >- 

>-  >* 

>- 

>-  >- 

4-4^4-4  4-4-. 

44  C  U-> 

00  : 

54 

33  S)  «-  1/ 

<0  >M  <» 

'OH!  O 

o  o 

O  n 

o  n 

o  o 

in  in 

«n 

©  © 

1  U  M  V  X 

1  v  V  w  <0 

rs 

45 

n  x 

(•■s  n 
X  x 

rs  n 
X  X 

rs  os 
x  x 

n  n 
x  x 

in 

X 

os  rs 

X  X 

1  4  4J  <0  V  OS 

l  v  <0  y  < 

l  U  u  ^4  <0  zo 

— 1  X  44  < 

j  c  x  c  \ 

e  e  —  a;  «w 

4-1 

V 

X 

n 

CM 

04 

<n 

04 

y 

04 

CO 

y 

04 

y 

rs 

y 

3  —  3  ^ 

,  9  0-* 

1  o  O  ©  <0  3) 

!  U  ©  <0S  >  — 

a 

O '  -* 

CSl 

04  04 

CO  04 

-y  cm 

m  cm 

X  04 

rs.  in 

oo  n 

!  33  II  ft  33  3) 

u  . 

V 

a. 

z, 

1 

•30 

r-s 

X 

X 

X 

X 

X 

y 

X 

y 

X 

lx-  -  X  X 
;  H  M  M  f-  H 

X 

s/ 

w 

w 

w 

w 

1  (0  XU 

Cti 


23 


t 

\ 

< 

s-V- 


I  tUlC tflia  ‘faVftrifa  rii&rftJ'miUtf  ti  Wi'aliU 


v,  .* 


W.' 


*■, 


-■  .'■'  !'  -■  >  ,%  v  ^\-  ^v"v“C*“3 


The  back  surface  gage  configuration  shown  in  Figure  2  was  used  in 
all  the  experiments.  A  flyer  plate  inclined  at  an  angle  of  90°  -  9  to 
the  direction  of  projectile  motion  impacts  a  parallel  specimen  plate 
backed  by  the  LiNbO^  gage.  To  distinguish  from  the  laboratory 

I  I 

coordinate  system  (X^-system),  we  define  a  X^-system  with  the  X2~axis 

t 

parallel  to  the  gage  thickness  direction,  and  the  X^-axis  coincident 
with  the  X^-axis. 

For  a  general  impact  (9  *  0°),  two  waves  are  propagated  into  the 
flyer  plate  and  the  sample  plate:  a  compressive  wave  and  a  shear 
wave.  The  compressive  wave,  which  travels  at  a  faster  velocity  through 
the  sample  plate,  impacts  the  specimen-gage  interface,  producing  a 
transmitted  wave  in  the  LiNbO^  and  a  reflected  wave  in  the  sample 
plate.  The  slower  traveling  shear  wave  in  the  sample  plate  interacts* 
with  the  reflected  compressive  wave  before  arriving  at  the  specimen-gage 
interface . 

The  use  of  a  back  surface  gage  for  determining  gage  calibrations  is 
not  an  optimal  arrangement  because  of  the  need  for  an  accurate  knowledge 
of  the  shear  response  of  the  specimen,  the  wave  interactions  in  the 
specimen,  and  the  shear  strength  of  the  specimen-gage  interface.  The 
advantage  of  using  the  specimen  (or  buffer)  plate  is  that  it  provides  a 
temporal  separation  between  the  compression  and  the  shear  wave.  This 
separation  permits  an  evaluation  of  the  gage  response  to  each  wave  and 
is  useful  in  modeling  the  response  of  a  rotated-cut  crystal.  If  the 
intent  is  to  measure  the  piezoelectric  constants,  then  compression-only 
(9  -  0°)  experiments  with  no  buffer  plates,  such  as  those  performed  by 
Graham, ^  are  desirable.  In  our  experiments,  the  flyer  and  specimen 
plates  were  made  of  PMMA  (Rohm  and  Haas  Type  II  UVA  Plexiglass).  This 
material  has  been  studied  both  under  compression^  and  shear  loading. ^ 
The  specimen  plate  thickness  was  typically  1  mm. 


In  the  present  case,  we  assume  linearity  and  ignore  this  interaction. 


For  experiments  where  a  shear  wave  was  introduced  in  the  sample 
plate  (9  *  0°),  the  direction  of  the  shear  motion  was  always  along  the 

f 

Xj-axis.  The  two  coordinate  systems  shown  in  Figure  2  are  related  to 
the  crystal  orientation  in  the  following  manner.  The  rotated  Y-cut 
crystals  involve  rotations  about  the  crystallographic  X-axis.  We  denote 
a  new  set  of  axes  relative  to  the  disk  of  the  cut  crystal  by  X',  Y',  and 
Z'.  The  Y'-axis  is  normal  to  the  crystal  disk,  and  is  related  to  the 
crystallographic  Y-axis  by  a  rotation  of  163°  or  165.5°  about  the 
crystallographic  X-axis.  As  a  result  the  X'-axis  coincides  with  the 
crystallographic  X-axis,  with  the  Z'-axis  normal  to  it  in  the  plane  of 
the  disk.  The  supplier  had  marked  both  the  X'-  and  Z'-axes  on  the  sides 
of  the  crystals.  As  indicated  earlier,  the  optimal  gage  is  one  that 
produces  an  electrical  output  only  when  sheared  along  the  Z'-direction. 

In  Table  1,  Z'  =  0°  means  that  the  Z'-axis  marked  on  the  gage 

1 

housing  was  oriented  along  the  X^-axis;  Z'  »  90°  indicates  that  it  was 
aligned  with  the  X^-  and  X^-axes.  The  purpose  of  the  Z'  =  90°  alignment 
in  combined  compression  and  shear  experiments  was  to  check  the  cross¬ 
axis  shear  effects.  For  this  alignment,  the  X'-  (or  X-)  axis  is  along 

» 

the  X-j-axis,  which  is  the  direction  of  shear.  A  signal  for  this 

alignment  would  imply  a  cross-axis  shear  effect.  Most  of  our 

experiments  were  compression-only  (9  =  0°)  experiments.  Aligning  the 

»  * 

Z'-axis  along  either  the  X3-  or  X^-axis,  was  advantageous  because  it 
permitted  us  to  evaluate  the  effect  of  impact  tilt,  as  discussed  later. 

The  impact  experiments  were  performed  using  the  SRI  gas-gun 

facility  for  studying  large-amplitude  one-dimensional  compression  and 

26 

shear  wave  propagation  in  solids.  Details  of  the  gun  barrel,  target 
holder,  and  related  instrumentation  for  measuring  particle  velocities, 
projectile  velocities,  and  impact  misalignment  are  presented  in 
Reference  26. 

The  target  assembly  for  the  impact  experiments  was  constructed  by 
bonding  the  LiNbO-j  gage  to  a  PMMA  specimen  plate  that  was  typically  5  mm 
thick.  Epon  R15  epoxy  resin  and  hardner  were  used  for  bonding,  and  the 
epoxy  layer  was  made  very  thin  (less  than  1-2  pm).  The  PMMA  specimen- 


gage  assembly  was  potted  in  aluminum  target  ring  as  shown  in  Figure  3. 
The  front  face  of  the  target  assembly  was  machined  and  ground  to  achieve 
a  1  mm  thickness  for  the  PMMA  plate.  The  completed  projectile  and 
target  assembly  shown  in  Figure  3  are  for  a  single-gage,  compression- 
only  experiment. 

The  U-shaped  vapor-deposited  gage  on  the  PMMA  surface  is  used  to 
measure  the  particle  velocity  at  the  impact  surface.  For  compression- 
only  experiments,  this  measurement  is  not  necessary  because  the  particle 
velocity  is  one-half  the  projectile  velocity.  However,  in  combined 
compression  and  shear  experiments,  this  measurement  is  essential  for 
determining  the  shear  particle  velocity  at  the  impact  face.  The  three 
sets  of  vapor-deposited  tilt  pins,  at  120°  intervals,  are  used  to 
measure  the  impact  misalignment.  The  fourth  set  of  metallic  leads  are 
0.001-inch-thick  copper  foils  used  to  trigger  the  electronics.  The  long 
outrigger  on  the  projectile  is  necessary  for  the  particle  velocity 
measurements,  as  indicated  in  Reference  26.  The  target  assembly  for  the 
combined  compression  and  shear  experiments  reported  here  is  identical  to 
that  for  the  compression-only  experiments.  The  projectile  assembly  is 
also  similar  except  that  the  flyer  plate  is  mounted  by  cutting  the 
outrigger  at  an  angle.  All  components  are  precisely  aligned  using  the 
three  holes  in  the  target  ring,  shown  in  Figure  3. 

When  we  first  began  these  experiments,  we  measured  the  tilt  only  to 
ensure  that  the  impact  misalignment  was  small.  However,  we  later  found 
that  the  tilt  measurements  were  of  particular  importance  because  an 
accurate  quantitative  analysis  of  the  data  needs  to  incorporate  the 
effects  of  tilt,  as  described  in  Section  4.3.  We  recommend  that  in 
future  experiments  the  impact  misalignment  be  measured  with  an  even 
higher  degree  of  accuracy. 
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Figure  3.  The  projectile  and  target  assembly  used  in  the  impact  experiments. 

(The  vapor-deposited  metallic  coating  on  the  projectile  shorts  the 
vapor-deposited  tilt  pins  on  the  target  to  provide  a  measurement 
of  the  impact  alignment.) 
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3.2  EXPERIMENTAL  RESULTS 


Table  2  summarizes  the  eight  experiments  performed,  which  involved 
14  gages.*  Of  these  experiments,  only  Experiment  1  is  taken  from  our 
feasibility  study  reported  earlier.  The  voltage-time  profiles  from  all 
eight  experiments  are  shown  in  Figures  4  and  5  and  are  discussed 
individually  in  this  subsection.  Appendix  F  shows  more  detailed 
current-  or  voltage-time  plots  prepared  after  digitizing  the  data. 

In  Table  2  the  impact  angle  refers  to  the  value  of  9  as  defined  in 
Figure  2(a).  Impact  tilt  is  based  on  the  closure  time  across  the  three 
sets  of  pins  shown  in  Figure  3.  Although  the  absolute  values  of  the 
tilt  measurements  are  probably  accurate  to  within  only  10  to  15  percent, 
the  values  cited  in  Table  2  are  good  representations  of  the  relative 
differences  between  the  shots.  The  fifth  column  of  Table  2  lists  the 
orientation  of  the  normal  to  the  line  of  contact  during  closure  relative 
to  the  X^-axis  as  defined  in  Figure  6.  For  example,  a  nonzero  value  of 
the  tilt  in  Experiment  3  implies  that  the  impact  was  not  simultaneous 
over  the  entire  face  of  the  specimen  plate;  the  90°  value  in  the  fifth 
column  indicates  that  the  closure  between  the  flyer  and  sample  plate 
proceeded  along  the  X^-axis.  Effectively,  this  results  in  a  shear 
deformation  along  the  X^-direction  even  for  a  compression-only  (9  *  0°) 
experiment.  As  discussed  later,  this  effect  is  significant  and  needs  to 
be  accounted  for  in  the  data  analysis. 

The  next  three  columns  of  Table  2  are  defined  as  in  Table  1.  The 
measured  current,  in  column  nine,  is  calculated  using  the  voltage  value 
at  the  first  jump  and  the  effective  resistance  values  cited  in  Table  1. 
In  using  piezoelectric  gages,  it  is  generally  preferable  to  use  current 
calibration  steps  rather  than  voltage  calibration  steps;  this  procedure 


Several  other  experiments  were  performed,  but  they  did  not  yield 
meaningful  data,  and  the  results  are  not  listed  here.  In  these 
experiments  either  the  gage  broke  or  foam  impactors  were  used.  The 
intent  of  using  the  foam  impactor  was  to  obtain  low  stresses  without 
lowering  the  impact  velocity. 
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avoids  measuring  Che  effective  resistance  separately.  We  believe  our 
calibration  procedure  is  reasonable  (within  1-2  percent)  because  the 
effective  resistance  was  carefully  measured.  In  fact,  the  largest  error 
occurs  in  reading  the  voltage  jump.  Because  of  the  significant  current 
ramping  following  the  initial  jump  and  because  of  the  rounding,  seen  in 
most  of  the  records,  it  is  difficult  to  identify  the  first  jump.  This 
error  is  significant  for  Experiment  6  as  seen  in  Figure  5(f).  The  three 
current  values  shown  in  parenthesis  in  the  current  column  correspond  to 
the  jump  following  the  shear  signal  and  include  the  contributions  from 
both  the  compression  and  shear  waves. 

The  duration  of  the  signal,  shown  in  Figures  4  and  5,  permits  a 
measurement  of  the  quasi-longitudinal  wave  (shock)  velocity  in  the 
LiNbOj  crystals,  which  is  indicated  in  the  second-to-last  column. 
Averaging  the  five  163°  Y-cut  values  and  the  three  165.5°  Y-cut  values, 
we  obtain  the  following  wave  velocities: 

163°  Y-cut:  6.73  mm/ ps 

165.5°  Y-cut:  6.72  mm/ ps 

Using  the  C^j  values  derived  from  Smith  and  Welsh's  paper, ^  we 
calculated  the  quasi-longitudinal  wave  velocity,  as  outlined  in  Appendix 
B.  The  calculated  wave  velocities  for  these  two  orientations  are  within 
0.1  percent  of  each  other,  and  we  cite  an  average  calculated  value  of 
6.743  mm/ ps  for  both  orientations.  Given  an  experimental  scatter  of  ±  1 
percent  in  the  data,  the  measured  wave  velocities  agree  with  the 
calculated  wave  velocities.  To  be  consistent  with  our  numerical 
calculations,  we  used  the  6.743  mm/us  value  in  analyzing  the  current 
data. 

In  Section  2.1,  the  current  output  for  the  idealized  situation  was 


derived  as 
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Figure  4.  Voltage-time  profiles  from  experiment  1  (78-2-46) 
with  163°  Y-cut  (Z'  =  0°)  LiNbOg  gage.  [The  top 
trace  represents  output  from  the  LiNb03  gage  (1  V/div) 
The  bottom  trace  represents  transverse  particle  velocity 
imparted  to  the  PMMA  (0.05  V/div).  Time  scale  is 
common  to  both  traces  (0.2  /ns/div) ] . 
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(c)  4(81-2-23) 
163°  (Z'  =  0°) 


(b)  3(81-2-22) 

UpperTrace:  165.5°  (Z' =  0°) 
Lower  Trace:  163°  (Z'  =  0°) 


(d)  4(81-2-23) 
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Figure  5.  Voltage-time  profiles  from  experiments  2  through  8 
with  rotated  Y-cut  LiNbO,  gages.  [The  voltage  and 
time  scale  per  division  are  shown  in  the  photographs. 
The  voltage  scale  in  (a)  was  3  V/div.] 


(e)  5  (81-2-24) 
163°  (Z‘  =  0°) 


(f)  6(81-2-431 
165.5°  (Z'  =  0°) 


<g)  6  (81-2-43) 
165.5°  (Z'  =  90°) 


Figure  5.  Voltage-time  profiles  from  experiments  2  through  8 
with  rotated  Y-cut  LiNb03  gages.  [The  voltage  and 
time  scale  per  division  are  shown  in  the  photographs. 
The  voltage  scale  in  (e)  was  3  V/div.  The  voltage  and 
time  setting  in  (f)  and  (g)  are  alike.]  (Continued). 
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(i)  8  (84-543) 

Upper  Trace:  163°(Z'  =  0°) 
Lower  Trace:  163°(Z'  =  90°) 


Figure  5.  Voltage-time  profiles  from  experiments  2  through  8 
with  rotated  Y-cut  LiNb03  gages.  (The  voltage  and 
time  scale  per  division  are  shown  m  the  photographs.) 
(Concluded). 


Figure  6.  Impact  closure  as  seen  from  the  projectile  side. 

[The  X[  axes  are  the  same  as  shown  in  Figure  2(b). 
The  points  1 ,  2,  and  3  refer  to  the  tilt  pins.  The 
two  axes  labeled  Z'  represent  the  orientations  of 
the  projection  of  the  Z'-axis  onto  the  plane  of  the 
crystal  disk  in  the  two  gage  orientations  typically 
used  in  our  experiments.  For  a  nonsimultaneous 
closure,  the  line  of  contact  is  defined  by  the  angle  <t>.) 
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where  i  is  the  gage  thickness,  A  is  the  eftective  electrode  area,*  Uc  is 
the  longitudinal  wave  velocity,  and  Ug  is  the  shear  wave  velocity. 

Using  the  data  given  in  Table  1  and  the  value  of  Uc  cited  above,  we 
calculated  the  compressive  polarization  (P  ),  listed  in  the  last  column 
of  Table  2.  The  polarization  due  to  shear  (P  )  is  shown  in  parenthesis 
for  the  three  compression-and-shear  experiments  along  the  163°  Y-cut 
orientation.  A  calculated  Ug  value  of  4.58  mm/|is  was  used  in 
determining  the  Pg  value.  Note  that,  unlike  the  current  values  in 
parenthesis,  the  shear  polarization  values  are  due  only  to  the  shear 
deformation.^  That  is,  the  shear  polarization  values  were  calculated 
from  the  difference  in  the  current  values. 

The  method  used  to  analyze  the  current  data  at  the  initial  jump  in 
our  experiments  is  similar  to  that  used  by  Graham  for  his  LiNbO^ 
data. ^  The  one  minor  difference  is  that  Graham  incorporates  finite 
strain  correction  in  his  results.  In  our  experiments,  the  strain  values 
are  less  than  0.2  percent,  and  this  correction  is  ignored  in  calculating 
the  polarization. 

As  mentioned  earlier,  the  effect  of  tilt  is  an  important  factor. 

As  shown  in  Figure  6,  not  only  is  the  tilt  angle  (column  4  in  Table  2) 
important,  but  the  orientation  of  the  line  of  closure  with  respect  to 

I 

the  X^-axis  (column  5  in  Table  2)  also  influences  the  result  in  two 
ways.  First,  the  angle  <t>  determines  the  relative  magnitude  of  the  shear 
strains  (  =  an<*  (  =  ^[3);  second,  for  two-gage  experiments, 

the  time  correlation  between  the  two  sets  of  signals  depends  on  the 
angle  A  general  theoretical  analysis  for  quantitatively 

incorporating  the  effects  of  impact  tilt  by  assuming  the  LiNbO^  as  an 
isotropic,  elastic  plate  is  presented  in  Appendix  E.  Further  discussion 
of  tilt  effects  using  this  analysis  is  presented  in  Section  4. 


*A  is  computed  using  a  diameter  that  extends  to  the  center  of  the 
insulating  gap  in  conformity  with  Graham’s  method. ^ 

^This  is  only  approximately  true  because  we  have  not  subtracted  the 
small  ramping  contribution  of  the  compression  wave. 


Results  from  the  Individual  experiments  (Figures  4  and  5)  are 
discussed  qualitatively  below.  These  data  are  analyzed  and  the 
piezoelectric  constants  are  calculated  in  Section  4. 

2 

Experiment  1  (78-2-46):  Figure  4,  taken  from  our  earlier  work, 
shows  the  shear  particle  velocity  at  the  impact  surface  along  with  the 
LiNbO^  output  on  a  common  time  scale.  The  measured  shear  particle 
velocity  was  approximately  90  percent  of  the  value  expected  from  the  "no 
slip"  condition.  The  LiNbi^  signal  shows  a  small  signal  from  the 
compression  wave  and  a  larger  signal  for  the  shear  wave.  We  did  not 
have  the  precise  resistance  values  and  hence  the  current  values  cited  in 
Table  2  are  good  to  only  4  or  5  percent.  However,  the  ratio  of  the 
compression  signal  to  shear  signal  does  not  depend  on  the  resistance  and 
should  be  fairly  accurate.  On  the  basis  of  the  tilt  angle  and 
the  <t>  value  in  Table  2,  we  expect  some  shear  contribution  even  for  the 
compression  wave. 

Experiment  2  (81-2-15):  This  experiment  was  designed  to  address  an 
important  objective  of  this  study:  the  effect  of  cross-axis  shear.  The 
two  gages  were  oriented  with  their  Z'-axis  parallel  (Z'  =  0°)  and  per¬ 
pendicular  ( Z *  ■  90°)  to  the  X3~axis.  In  Figure  5(a),  a  clear  two-step 
signal,  similar  to  that  in  Figure  4,  is  seen  for  the  Z'  *  0°  gage;  the 
Z'  =  90°  gage  does  h-'*-  show  a  two-step  signal.  This  result  confirms  the 
theoretical  predictions  that  the  shear  gage  should  display  no  cross-axis 
effect.  Because  of  the  direction  of  the  impact  closure  line  (<t>  *  140°), 
the  Z'  *  0°  gage  signal  began  earlier,  as  expected.  There  are  two 
features  of  this  experiment  that  we  believe  are  caused  by  the  impact 
tilt:  the  difference  in  the  compression  signal  between  the  Z'  =  90°  and 

the  Z'  *  0°  gage  (see  Table  2),  and  the  small  decrease  in  the  Z'  *  90° 
signal  at  the  time  of  the  shear  wave  arrival  at  the  specimen-gage  inter¬ 
face.  We  address  these  issues  in  the  next  section. 

Experiment  3  (81-2-22):  This  experiment  was  designed  to  compare 
the  output  of  163°  and  165. 5 0  Y-cut  crystals  in  a  compression-only 
experiment.  Such  an  experiment  is  also  useful  in  checking  the  piezo¬ 
electric  constants  because  the  165.5°  Y-cut  is  a  specific  direction 
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on  the  basis  of  the  literature  values.  Fortunately,  this  experiment 
also  had  the  lowest  tilt  and  <t>  was  90°.  Thus,  the  gages  were  subjected 
to  identical  loading.  In  Figure  5(b),  the  initial  current  output  of  the 
165.5°  Y-cut  gage  is  higher  than  that  of  the  163°  Y-cut  gage  (also  see 
Figure  F.3).  As  expected  for  =  90°,  both  gage  signals  begin  at  nearly 
the  same  time.  The  165.5°  Y-cut  gage  signal  actually  begins  earlier 
even  though  it  had  a  thicker  specimen  plate.  This  suggests  that  the 
correlation  with  our  tilt  measurements  is  reasonably  good  but  not 
exact.  For  an  exact  correlation,  the  centers  of  the  two  gages  should  be 
emplaced  comp' otely  symmetric  to  the  tilt  pins. 

The  signals  drop  at  a  time  that  corresponds  to  the  arrival  of  a 
shear  wave  at  the  specimen-gage  interface.  Because  of  the  higher 
voltage  sensitivity  for  the  163°  Y-cut  gage,  this  drop  is  more 
discernible  in  the  record.  Because  no  shear  deformation  was  introduced 
into  the  specimen,  the  small  shear  wave  generated  due  to  tilt  (see 
Appendix  E)  may  be  the  cause  of  this  drop  in  the  signal.  An  estimate  of 
this  drop  is  made  in  the  next  section.  (This  is  the  same  as  the  tilt 
correction) . 

Experiment  4  (81-2-23):  This  experiment  was  performed  to  confirm 
the  findings  of  Experiment  3  and  provide  data  at  a  different  compressive 
level  for  calculating  the  piezoelectric  constants.  Unfortunately,  the 
impact  tilt  was  nearly  2.5  times  the  tilt  in  Experiment  3  although  <t>  was 
90°.  Because  <t>  was  90°,  we  expected  the  two  gage  signals  to  begin 
simultaneously.  The  signal  of  the  163°  Y-cut  gage  is  observed  to  start 
slightly  early,  confirming  our  remarks  earlier  concerning  the  limited 
accuracy  of  our  tilt  measurements.  The  overall  features  of  the  gage 
records  are  similar  to  those  shown  in  Figure  5(b)  for  Experiment  3.  The 
drop  in  the  signal,  discussed  above,  occurs  in  both  gage  records. 
Linearly  scaling  the  results  of  Experiment  3,  we  expected  a  current  jump 
of  82  mA  for  the  165.5°  Y-cut  gage.  Instead,  the  observed  jump  was  35 
mA.  The  current  jump  for  the  163°  Y-cut  gage  was  5  to  7  percent  higher 
than  the  jump  in  Experiment  3  although  linear  scaling  suggests  a 
decrease  of  30  percent.  Again,  the  influence  of  tilt  is  at  least  a 
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qualitative  explanation  for  this  because  it  would  have  a  larger  relative 
effect  for  the  163°  Y-cut  orientation.  However,  without  a  quantative 
analysis  this  is  only  conjecture. 

Experiment  5  (81-2-24):  This  experiment  was  designed  to  provide 
data  similar  to  that  from  Experiment  2  but  at  a  higher  stress. 
Unfortunately,  the  tilt  was  40%  higher  than  in  Experiment  2,  and  the  Z’ 

=  90°  gage  did  not  produce  any  meaningful  data.  The  Z'  =  0°  gage  again 
shows  the  characteristic  two-step  signal  expected  for  combined 

compression  and  shear  loading.  The  shear  polarization  value,  Ps»  scales 
very  well  between  this  experiment  and  Experiment  2.  This  value, 
measured  as  the  difference  between  the  two  current  steps,  is  largely 
unaffected  by  tilt.  The  current  jump  due  to  compression  for  this 
experiment  is  larger  than  expected  from  linear  scaling  with  Experiment 
2.  However,  the  tilt  value  in  Experiment  5  was  higher  than  that  in 
Experiment  2  (the  <J>  value  was  comparable)  and  this  can  influence  the 
results . 

Experiment  6  (81-2-43):  This  experiment  was  designed  to  examine 
the  differences  in  the  effect  of  tilt,  if  any,  in  a  compression-only 
experiment.  Both  gages  had  the  165.5°  Y-cut  orientation.  One  was 
aligned  with  Z'  *  0°  and  the  other  with  Z'  =*  90®.  As  seen  in  Figures 
5(f)  and  5(g),  the  two  voltage-time  profiles  are  quite  different.  We 
have  spent  considerable  time  in  trying  to  understand  this  result  but 
have  had  little  success.  The  marked  roundedness  of  the  initial  jump  in 
both  records  adds  uncertainty  to  a  comparison  of  their  magnitudes,  and 
the  profiles  beyond  the  first  jump  are  markedly  different.  One  is 
forced  to  conclude  that  tilt  and  variations  in  gages  must  be  the  cause 
because  there  is  no  other  way  to  distinguish  between  these  gages  in  a 
compression-only  experiment.  The  impact  tilt  and  the  <J>  value  in  Table  2 
suggest  that  the  shear  generated  by  the  tilt  will  increase  the 
polarization  for  the  %'  =  0°  gage.  The  Z'  *  0°  gage  profile  from  this 
experiment  is  comparable  to  the  165.5°  Y-cut  (Z'  *  0°)  gage  profile  from 
Experiment  3,  shown  in  Figure  5(b).  The  higher  initial  jump  in 
Experiment  6  compared  with  the  initial  jump  in  Experiment  3  also 


suggests  tilt  effects.  It  is  the  Z'  =  90°  gage  profile  in  this  experi¬ 
ment  that  remains  a  mystery. 

•fa 

Experiment  7  (84-542)  :  This  experiment  was  performed  to  reproduce  and, 
we  hoped,  provide  understanding  of  the  response  observed  for  the  165.5° 
(Z'  =  90°)  gage  in  Experiment  6.  The  measured  voltage-profile  shown  in 
Figure  5(h)  is  considerably  different  from  the  profile  in  Figure  5(g). 
Because  the  tilt  for  this  experiment  was  small,  the  rise  time  in  Figure 
5(h)  is  smaller  than  Experiment  6.  The  tilt  result  from  this  experiment 
is  comparable  to  Experiment  3  and  the  overall  profiles  in  Figure  5(d) 
and  Figure  5(g)  are  similar.  The  measured  current  jump  of  92.5  mA  is 
markedly  different  from  the  values  in  Experiment  6  but  is  closer  to  the 
current  values  in  Experiments  3  and  4.  Before  drawing  any  final  conclu¬ 
sions,  the  tilt  analysis  needs  to  be  incorporated.  On  the  basis  of  this 
experiment,  and  the  similarity  with  Experiments  3  and  4,  we  conclude 
that  the  results  from  the  165.5°  (Z'  =  90°)  gage  in  Experiment  6  are 
anomalous.  The  cause  of  the  anomoly  is  presently  unknown. 

Experiment  8  (84-543)*:  Because  the  results  from  Experiment  6  were  not 
useful,  we  performed  a  compression-only  experiment  with  Z'  =  0°  and 
Z'  =  90°  gages  to  examine  the  effects  of  tilt.  In  the  absence  of  tilt 
these  gages  should  provide  identical  response.  The  only  gages  that  were 
available  to  us  from  the  original  batch  were  163°  Y-cut  crystals.  The 
measured  voltage  profiles  are  shown  in  Figure  5(i).  As  indicated  in 
Table  2,  the  current  and  polarization  at  the  first  jump  are  50  percent 
different.  This  difference  can  be  due  to  two  factors:  tilt  effects  or 
variations  in  gage  orientations.  These  dat  \  are  analyzed  in  the  next 
section  and  compared  with  theoretical  predictions  incorporating  tilt 
effects  to  determine  the  cause  of  the  observed  difference. 


Experiments  7  and  8  were  performed  at  the  Shock  Dynamics  Laboratory 
of  Washington  State  University.  The  assistance  of  Pete  Majewski  and 
Martin  Williams  is  gratefully  acknowledged . 


The  experimental  results  from  all  the  experiments  are  summarized  as 
follows : 

(1)  The  163°  Y-cut  LiNbO^  crystals  show  considerably  larger 
sensitivity  to  shear  loading  than  to  compressive  loading. 

(2)  The  163°  Y-cut  LiNbC^  does  not  appear  to  be  sensitive  to 
cross-axis  shear. 

(3)  There  is  considerable  experimental  evidence  to  suggest  that, 
because  of  the  large  shear  sensitivity  of  these  gages,  we  need 
to  account  for  tilt  effects  in  data  analysis.  Not  only  is  the 
magnitude  of  the  tilt  angle  important,  but  the  relative 
orientation  of  the  impact  closure  line  can  also  make  a 
difference. 

(4)  The  response  of  the  crystal  in  the  Y-cut  165.5°  orientation  to 
compression  loading  is  considerably  higher  than  the  response 
of  the  crystal  in  the  163°  Y-cut  orientation. 

(5)  The  results  of  Experiment  6,  in  which  two  165.5°  Y-cut  gages 
(Z*  =0°  and  Z'  *  90°)  were  subjected  to  compression-only 
loading,  are  not  currently  understood.  Based  on  the  results 
of  Experiment  7,  the  Z*  *  90°  result  in  Experiment  6  is 
believed  to  be  anomalous. 

(6)  The  large  current  rise  beyond  the  initial  jump  suggests  the 
importance  of  modeling  electromechanical  coupling  at  late 
times . 

(7)  Some  of  the  experimental  results  do  not  scale  with  projectile 
velocity,  suggesting  effects  of  tilt  or  variability  of  the 
orientation  of  the  crystals.  The  latter  effect  would  be 
particularly  important  for  the  163°  orientation. 

(8)  The  wave  velocity  measurements,  inferred  from  the  duration  of 
the  signal,  show  ail  percent  scatter.  Within  this  scatter, 
these  values  match  the  calculations. 


SECTION  4 


ANALYSIS  AND  DISCUSSION 

This  section  compares  the  experimental  results  with  the  theoretical 
predictions  using  the  uncoupled  elastic-dielectric  model  outlined  in 
Sections  2.1  and  2.2.  We  believe  that  the  uncoupled  approximation  is 
appropriate  for  analyzing  the  initial  jump  in  our  experiments. 

Potential  causes  for  discrepancy  between  the  theory  and  experiment  are 
examined.  Tilt  effects,  using  the  theoretical  developments  presented  in 
Appendix  E,  are  estimated.  The  overall  findings  are  discussed  at  the 
end  of  the  section. 


4.1  NUMERICAL  SIMULATION  OF  THE  EXPERIMENTS  (MECHANICAL  VARIABLES) 


One-dimensional  wave  propagation  code  calculations  were  performed 
using  the  subroutine  ANELAS  described  in  Appendix  B  for  computing  the 
stress  in  anisotropic  LiNb03  crystals.  The  pertinent  mechanical 
variables  in  the  LiNbO-j  for  each  of  our  experiments  are  listed  in  Table 
3.  The  LiNbO-j  elastic  constants,  C?.,  used  in  these  calculations  are 

C  it  J 

from  the  work  of  Smith  and  Welsh.  The  PMMA  was  modeled  as  a  nonlinear 

elastic  material  using  the  following  material  properties: 


1 

Initial  density:  Pq“  1.185  g/cm 

2  3 

Pressure-volume  relation:  P  (kbar)  *  62  p  +  295  p  +  134  p 


Shear  modulus: 


where  p 
G  (kbar) 


»  - - 1 

Po 

=■  22.5  +  P 


(4-1) 


*  E 

These  authors  report  only  .  constants  in  their  paper.  These  con¬ 
stants  were  converted  to  CP.J using  the  procedure  indicated  in 
Appendix  8. 


Table  J.  Numerical  simulation  of  LiNbOj  experiments  using  C  *.  from  reference  S. 
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All  gages  listed  here  are  oriented  with  Z'  -  0°.  See  Table  1,  footnote  a. 
h'lle  results  for  these  experiments  can  be  scaled  from  the  results  for  Experiment  3  (81-2-22). 
*  The  two  sets  of  results  for  the  compression  and  shear  experiments  are  explained  In  the  text. 
^ The  figures  In  parentheses  refer  to  angle  defined  in  Figure  2. 


For  compressive  stresses  below  7.6  kbar,  the  results  of  this  model  are 

in  good  agreement  with  the  PMMA  data  of  Barker  and  Hollenbach, ^  and 
25 

Gupta.  In  all  the  numerical  simulations  the  PMMA  buffer  plate  was 
modeled  as  1  mm  thick,  which  is  slightly  different  from  the  actual 
thicknesses  listed  in  Table  2.  Because  of  the  elastic  nature  of  the 
calculation,  this  difference  does  not  affect  the  wave  amplitudes. 
Comparison  with  arrival  times  will  require  appropriate  scaling.  The 
results  shown  in  Table  3  are  based  on  zero  tilt,  that  is,  perfect 
alignment  at  impact.  Hence  tilt  effects  need  to  be  separately  accounted 
for. 

For  compression-only  experiments  in  Table  3  there  is  only  a  single 
wave  in  the  PMMA  buffer  plate.  The  stresses  and  strains  are  listed  in 
the  matrix  notation.  The  results  scale  approximately  with  the 
projectile  velocity;  the  scaling  is  not  exact  because  of  the  nonlinear 
response  of  the  PMMA.  For  the  165.5°  Y-cut  orientation,  the  LiNb02 
strain  response  is  nearly  uniaxial.  For  the  163°  Y-cut  orientation,  the 
shear  strain,  although  small,  can  have  significant  contributions  to  the 
polarization,  as  discussed  later.  Because  of  the  closeness  of  impact 
velocities,  the  results  for  Experiments  6  through  8  can  be  scaled  from 
the  results  for  Experiment  3. 

The  calculated  results  for  the  combined  compression  and  shear 
experiments  in  Table  3  were  results  using  a  no-slip  condition  at  both 
the  impact  and  the  PMMA/LiNbO^  interface.  The  two  sets  of  values  for 
each  experiment  result  from  the  two  waves  generated  in  the  PMMA  buffer 
plate.  The  first  entry  for  each  experiment  represents  the  LiNbO^  state 
after  the  compression  wave  arrival  at  the  PMMA/LiNbO^  interface;  the 
second  entry  represents  the  final  state  after  the  arrival  of  the  shear 
wave  at  the  gage  interface. 

The  shear  strain  e’  and  the  corresponding  stress  and  particle 
o 

velocity  are  zero  in  Table  3.  This  result  differs  slightly  from  the 
numerical  results  presented  in  Table  C-2  of  Appendix  C  for  Y-cut  and 
36°-rotated  Y-cut  crystals.  The  cause  for  this  slight  difference  is  the 
use  of  an  X-cut  quartz  crystal  as  impactor  in  contrast  to  the  use  of 
PMMA  as  an  impactor. 
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In  Appendix  B,  we  present  the  C?.  values  from  the  work  of  Smith  and 

5  4  ^  D 

Welsh  and  Warner  et  al.  Except  for  the  value  of  values  for 

the  other  constants  are  quite  close.  To  determine  the  effects  of 
scatter  in  the  C^j  values  on  our  results,  we  performed  a  series  of 
calculations  for  compression-only  experiments  using  the  values  of 
Warner  et  al.^  The  results  are  presented  in  Table  4.  Comparing  Tables 
3  and  4  shows  that  the  values  of  the  compressive  variables  are  almost 
identical.  However,  the  values  of  the  shear  variables  are  different, 
and  the  larger  relative  differences  observed  for  the  165.5°  Y-cut 
orientation  indicate  that  the  two  sets  of  constants  give  rise  to  some¬ 
what  different  specific  directions.  However,  both  sets  of  constants 
show  that  the  strain  response  along  the  165.5°  Y-cut  is  very  nearly 

uniaxial.  The  e!  differences  for  163°  Y-cut  orientation  are  not 
4 

negligible  and  provide  a  measure  of  one  source  of  scatter  in  comparing 
polarization  values  between  theory  and  experiment.  This  is  discussed  in 
subsection  4.3. 

4.2  ANALYSIS  OF  THE  POLARIZATION  MEASUREMENTS 

The  polarization  along  the  gage  thickness  for  a  rotated  Y-cut 
crystal  can  be  written  in  the  matrix  notation  as 

p2  “  e22  e2  +  e24  e4  (2-23) 

t  » 

There  is  no  contribution  from  the  e2^  term  even  when  is  nonzero 
because  for  a  Y-cut  or  a  rotated  Y-cut  LiNbO-j  crystal  is  zero. 

I  f 

Expressions  for  e22  and  e2^  are  presented  in  Appendix  D;  the  general 

f 

expression  for  P2  in  terms  of  e^j  is  shown  in  equation  (C.l). 

Using  the  strain  values  cited  in  Table  3,  the  piezoelectric 

I 

constants  e^j,  and  equation  (2-23),  we  can  compare  the  calculated  and 
measured  polarization  values,  as  shown  in  Table  5.  In  comparing  these 
values,  the  following  stipulations  are  important: 

(1)  The  measured  current  was  converted  to  polarization  using  the 
uncoupled  elastic-dielectric  analysis  of  Section  2.2  because 
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ignoring  the  electromechanical  coupling  is  reasonable  for  the 
initial  jump.  However,  we  have  had  to  also  ignore  the 
contribution  of  the  current  ig  in  equation  (2-24)  because  we 
have  no  way  of  assessing  it.  This  simplification  of  the 
analysis  of  Section  2.2  is  equivalent  to  neglecting  the  off- 
diaganol  elements  of  the  permittivity  tensor,  or  the  E  field 
components  normal  to  the  gage  thickness  direction  (see 
equation  2-22).  Without  a  means  of  estimating  ig,  it  is 
difficult  to  assign  error  bars  on  the  measured  polarization 
value. 

(2)  The  polarization  values  listed  in  parenthesis  in  Table  5  are 
the  shear  contributions.  In  calculating  these  values,  we 
assumed  that  the  impact  surface  and  the  PMMA/LiNbOj  interface 
bond  can  transmit  the  shear  stresses. 

(3)  The  two  sets  of  calculations  reflect  the  two  sets  of  e^. 
values  inferred  from  the  literature.  The  main  discrepancy  is 
in  the  constant,  as  discussed  in  Appendix  C.  The  two  set 
of  values  used  in  calculating  ?2  fr°®  equation  (2-23)  are  as 


follows: 

Orientation 

r  2 
e22  ) 

e24 

163°  Y-cut 

-8  x  10“3  (Old) 

-3.76  (Old) 

-3.49  x  10"1  (New) 

-3.86  (New) 

165.5°  Y-cut 

-3.58  x  10"1  (Old) 

-3.85  (Old) 

-7.12  x  10-1  (New) 

-3.94  (New) 

These  values  are  discussed  in  greater  detail  in  Appendix  D. 

Comparison  of  the  calculations  with  the  experimental  measurements 
leads  to  the  following  conclusions: 

(1)  The  observed  shear  contributions  to  polarization  in  Experiments 
2  and  5  closely  match  the  theoretical  values. 

(2)  The  results  of  Experiment  2  suggest  that  there  is  no  cross-axis 
sensitivity. 

(3)  The  agreement  between  the  theory  and  experiment  is  not  good  for 
compression  loading.  Most  of  the  compression  loading  data 
agree  better  with  the  old  set  of  constants.  However,  some 


* 

Experiment  1  is  less  reliable,  as  Indicated  in  Section  3. 


agree  with  the  new  set  of  constants;  no  systematic  trends  are 
obvious.  The  calculated  compression  response  for  163°  Y-cut 
crystals  also  depends  markedly  on  the  Cj?  value  chosen.  For 
example,  using  the  values  in  Table  4,  trie  calculated 
polarization  for  Experiment  3  in  Table  5  is  0.36  and  not  0.44. 

(4)  These  inconsistencies  and  the  inconsistencies  among  the  data 
themselves  (dicussed  in  Section  3)  suggest  that  contributions 
from  the  following  sources  of  errors  may  be  important:  (a) 
finite  tilt  values,  (b)  deviations  from  the  desired  crystal 
orientation,  (c)  uncertainties  in  the  knowledge  of  the  piezo¬ 
electric  constants  and  C,^  *  (d)  the  extra  current  term  in 
the  rotated-cut  derivation  discussed  in  Section  2.2,  and  (e) 
the  contribution  from  electromechanical  coupling  due  to  tilt, 
which  can  introduce  an  error  in  the  measured  current  jump  e.g. 
Experiment  6.  An  estimate  of  the  contributions  from  the  first 
two  sources  is  presented  in  the  remainder  of  this  section.  The 
contribution  from  the  last  two  sources  can  only  be  assessed 
after  further  theoretical  work. 


4.3  ASSESSMENT  OF  TILT  EFFECTS 

We  were  unable  to  develop  a  simple  method  to  estimate  tilt 
effects.  Therefore,  the  analysis  described  in  Appendix  E  was  developed. 
However,  this  analysis  is  not  exact  because  it  treates  the  LlNbO^  as  an 
isotropic  solid.  Despite  this  limitation,  we  believe  that  the  results 
of  our  tilt  analysis  provide  a  quantitative  assessment.  By  using  the 
compression  and  shear  wave  velocities  that  are  pertinent  to  our  crystal 

f  I 

orientations,  we  can  obtain  reasonable  values  of  e 2  and  £4  in  the  LiNbO^ 
for  use  in  equation  (2-23).  Also,  the  crystal  orientations  are  fairly 
close  to  the  specific  direction. 

As  indicated  in  Appendix  E,  a  finite  tilt  results  in  a  compression 
and  a  shear  wave  inclined  at  an  angle  to  the  impact  surface.  The  faster 
travelling  compression  wave  reaches  the  PMMA/LiNbO^  interface  and 
produces  a  compressive  stress  (o^)  and  a  shear  stress  (o^).  We 
calculated  the  contribution  to  polarization  from  the  shear  stress 
(actually,  shear  strain).  To  convert  the  stresses  into  strains,  we  used 
the  following  moduli  for  LiNbO^:  Longitudinal  modulus  ■  2114  kbar  and 
Shear  modulus  *  986  kbar.  These  values  are  equivalent  to  a  longitudinal 
wave  velocity  of  6.75  mm/ps  and  a  shear  wave  velocity  of  4.61  mra/ps. 


Both  values  are  very  close  to  the  wave  velocities  calculated  by  using 
the  C?j  values  and  equation  (B.6). 

Because  we  are  approximating  the  LiNbO-j  as  an  isotropic  crystal, 
the  other  stresses  will  not  be  matched  correctly.  In  the  present 
analysis  this  is  not  important  because  the  polarization  contribution  due 

t 

to  tilt  arises  from  the  €4  strain  caused  by  the  finite  tilt.  For  both 
orientations  (163°  and  165.5°  Y-cut),  has  the  same  value  for  a 

I 

particular  impact  experiment.  To  calculate  the  polarization,  PT,  we 

t  2 

used  a  &2U  value  of  -3.85  C/m  .  This  value  is  reasonable  because  the 
» 

value  for  the  two  orientations  using  either  set  of  constants  ranges 
from  -3.76  to  -3.94  C/m2. 

In  the  analysis  presented  in  Appendix  E,  we  also  calculated  the 
effect  of  the  slower  travelling  shear  wave,  generated  at  the  impact 
surface,  after  it  reaches  the  PMMA/LiNbC^  interface.  The  effect  of  this 
shear  wave  is  not  included  in  this  analysis  of  tilt  effects  because  we 
are  Interested  only  in  calculating  the  polarization  at  the  first  jump. 

The  results  incorporating  the  effects  of  tilt  are  summarized  in 

» 

Table  6.  The  compressive  strains,  €2*  were  derived  from  the  oblique 
plate  analysis  described  in  Appendix  E,  using  the  longitudinal  modulus 
value  cited  earlier.  Although  these  values  are  not  used  for 

I 

polarization  calculations,  they  can  be  compared  with  the  e 2  values  cited 

in  Table  3.  The  shear  strains  £4  were  calculated  from  the  shear  stress 

» 

due  to  finite  tilt.  The  tilt  contributions  to  the  polarization,  P^, 
were  calculated  from 

pj  =■  -3.85  €4  C/m2  (4-2) 

f 

We  interpret  the  tilt  contribution  as  follows:  the  PT  values 
represent  a  contribution  to  polarization  that  is  present  in  the 


The  eigenvalues  for  the  163°  and  165.5°  Y-cut  orientation  are  within 
0.2  percent  of  each  other,  and  we  use  the  average  calculated  values: 
6.743  mm/^s  (P-wave)  and  4.58  mm/^s  (S-wave). 


experiments  but  not  included  in  the  theory.  Hence,  the  P^,  values  in 

Table  6  should  be  added  to  the  calculated  polarization  values  shown  in 

Table  5.  We  used  the  polarization  values  in  Table  5  based  on  the  old 

*  * 

e^j  constants  to  add  to  the  values.  This  calculated  total  is 
compared  with  the  experimental  measurements  in  the  last  two  columns  of 
Table  6.  Except  for  Experiments  1  and  3,  the  addition  of  the  tilt  con¬ 
tribution  improves  the  agreement  between  theory  and  experiment.  It  is 
noteworthy  that  most  of  the  difference  in  the  output  of  the  two  gages  in 
Experiment  8  was  due  to  tilt  effects.  Without  tilt  effects,  the  two 
gages  ( Z *  -0°  and  Z'  =  90°)  would  appear  identical  as  expected.  This 
result  underscores  the  anomalous  response  of  Experiment  6.  Although 
differences  remain  between  theory  and  experiment,  the  results  are 
encouraging. 

In  Table  6,  we  did  not  calculate  the  tilt  correction  for  the 
contribution  of  the  shear  strain  to  the  polarization  values  shown  in 
parenthesis  in  Table  5.  First,  we  wanted  to  analyze  only  the  initial 
jump;  subsequent  results  are  strongly  influenced  by  the  electromechani¬ 
cal  coupling.  Second,  the  S-wave  due  to  finite  tilt  (not  due  to  the 
shear  deformation  deliberately  introduced  by  making  9  *  0  in  Figure  2) 
on  arrival  at  the  PMMA/LiNb03  interface  reduces  the  shear  stress  (or 
strain  to  a  small  value.  This  effect  is  completely  masked  by  the 
arrival  of  the  main  shear  wave.  Also,  this  result  suggests  that  the 
tilt  error  is  self-compensating  for  the  shear  polarization  part.  Note 
that  remarks  about  tilt  corrections  at  late  times  cannot  be  considered 
rigorous  because  of  electromechanical  coupling  effects. 

Comparisons  of  the  measured  and  calculated  polarization  values 
after  accounting  for  tilt  effects  in  an  approximate  manner  are 
summarized  as  follows: 


Using  the  other  set  would  further  increase  the  discrepancy  between 
theory  and  experiment. 


(1)  The  tilt  contributions  to  polarization  are  significant  because 
of  the  large  shear  sensitivity  of  the  gage. 

(2)  In  six  out  of  eight  measurements,  the  incorporation  of  tilt 
effects  reduces  the  discrepancy  between  theory  and  experiment. 

(3)  The  incorporation  of  tilt  effects  indicates  that  the  old  set 
of  piezoelectric  constants  provides  a  better  match  to  our 
data.  The  old  set  of  constants  are  presented  in  Appendix  D. 

The  analysis  presented  here  is  not  exact.  An  improved  analysis 
requires  that  the  anisotropic  response  of  LiNbC^  be  accounted  for  in 
calculating  tilt  effects.  We  also  believe  that  our  tilt  measurements 
can  be  improved  in  future  experiments. 

4.4  DISCUSSION 

The  analyses  presented  in  this  section  have  answered  some  questions 
but  also  raised  new  ones. 

On  the  positive  side,  the  present  work  has  confirmed  the  large 
sensitivity  of  the  163°  Y-cut  orientation  to  shear  loading  relative  to 
compression  loading.  The  one  experiment  designed  to  examine  sensitivity 
to  cross-axis  shear  suggests  that  the  gage  is  not  sensitive  to  cross¬ 
axis  shear.  As  predicted,  the  compression  sensitivity  of  the  165.5° 
Y-cut  orientation  is  considerably  higher  than  that  of  the  163°  Y-cut 
orientation.  A  fairly  complete  theoretical  framework  has  been  developed 
to  quantitatively  analyze  the  gage  response  to  well-defined  loading. 
Detailed  comparisons  of  theoretical  predictions  and  experimental 
measurements  have  been  presented.  For  shear  loading,  the  theory  and 
experiment  show  good  agreement  (within  3  percent). 

On  the  negative  side,  quantitative  agreement  between  theory  and 
experiment  is  not  good  for  compression  loading  (10  to  50  percent).  Some 
of  the  inconsistencies  among  the  experimental  results  have  not  been 
completely  resolved.  The  results  from  Experiment  6  are  believe  to  be 
anomalous  (see  Figure  5(f)  and  5(g)  in  Section  3).  Similarily,  the 
differences  between  the  first  jump  for  the  two  gages  in  Experiment  2 
(Table  5)  are  difficult  to  understand. 


One  is  tempted  to  speculate  that  neglecting  of  the  current  ig  term 
[see  equation  (2-24)]  in  converting  the  current  values  to  polarization 
may  be  a  major  source  of  error.  However,  this  error  is  probably 
systematic  because  the  experimental  conditions  are  fairly  similar.  The 
random  tilt  variations  can  only  explain  part  of  the  discrepancy.  Hence, 
we  are  forced  to  conclude  that  the  variability  in  the  orientations  of 
the  crystals  may  be  causing  some  of  the  observed  inconsistencies.  The 
supplier  assured  us  that  the  orientations  were  within  0.2°.  Although  we 
cannot  check  the  gages  that  have  been  impacted,  the  orientations  of  the 
remaining  gages  from  that  batch  were  independently  measured.*  These 
results  indicated  that  the  scatter  in  orientation  ranged  between  0°  and 
±  1.5°.  To  us  this  scatter  seems  high  because  the  supplier^  of  these 
gages  had  considerable  expertise  in  crystal  orientation.  Further 
experimental  work  is  needed  to  unequivocally  resolve  this  issue,  as 
indicated  in  Section  5. 

There  is  also  the  question  of  the  two  sets  of  piezoelectric 

constants.  The  main  difference  is  in  the  e.^2  constant.  The  value  cited 

5  2 

in  the  literature  (from  ultrasonic  experiments)  is  2.43  C/m  .  Our 

analysis  of  Graham's  shock  wave  data,^  described  in  Appendix  C,  yields  a 

value  of  2.82  C/m  .  Because  the  new  set  of  constants  (using  e22  =  2.82) 

are  based  on  impact  results,  one  is  tempted  to  believe  these  results. 

However,  our  present  data  are  better  matched  by  using  the  old  set  of 

constants  (&22  =  2.43).  Clearly,  the  value  of  e22>  because  of  its  large 

influence  on  the  results,  has  to  be  established  accurately.  There  are 

two  deficiencies  in  the  analysis  of  Graham's  data^.  First,  the  extent 

of  tilt  effects  is  not  known.  Second,  the  conversion  of  the  current 

jump  to  polarization  does  not  account  for  the  two-dimensional  effect  (ip 


.  ■'."I 


rtm 


We  are  grateful  to  Mr.  R.  A.  Graham  of  Sandia  National  Laboratories 
for  these  measurements.  The  standard  deviation  was  estimated  to  be 
±0.5°  in  these  measurements. 

^ Now  deceased. 


term)  discussed  in  Section  2.2.  When  this  latter  effect  is  not  properly 
modeled,  the  shock  response  of  rotated-cut  crystals  is  suspect.  Hence, 
we  recommend  using  the  old  set  of  constants  as  defined  in  Appendix  D. 

I 

An  added  complexity  in  interpretation  arises  because  the  new  &22 
2 

value  of  -0.349  C/m  along  the  163°  Y-cut  orientation  is  comparable  to 

t 

the  0-22  value  based  on  the  old  constants  but  for  a  crystal  misorienta- 
tion  of  ±  0.2°  (see  Appendix  D) .  There  is  also  the  question  of 
difference  in  the  values  reported  in  References  4  and  5.  As 

indicated  at  the  end  of  subsections  4.1  and  4.2,  this  difference  can 
markedly  influence  the  compressive  contribution  to  polarization  for  the 

t 

163°-rotated  cut  crystals.  If  the  e  value  used  in  Equation  (2.23)  is 

A  t 

taken  from  Table  4  instead  of  Table  3,  the  value  for  the  compressive 
response  of  the  163 “-rotated  cut  crystals  would  be  decreased  by 
approximately  15-20  percent. 

It  is  obvious  from  the  above  remarks  that  improved  reconciliation 
between  theoretical  predictions  and  experimental  measurement?  will 
require  further  experimental  and  theoretical  work,  and  an  improved  know¬ 
ledge  of  the  piezoelectric  constant  &22  and  ttie  Mastic  constant  C^. 
Because  of  the  many  potential  sources  of  errors,  it  is  hard  to  proceed 
further  without  additional  information.  The  need  for  resolving  the 
differences  between  the  theoretical  measurements  and  experimental 
results  for  developing  a  gage  package  and  the  specific  steps  to  resolve 
these  differences  are  discussed  in  the  next  section. 


SECTION  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  objectives  of  this  work,  outlined  in  Section  1.2,  have  largely 
been  met.  A  fairly  complete  theoretical  framework  has  been  developed  to 
analyze  the  response  of  rotated-cut  LiNbO^  gages.  The  experimental 
work,  designed  to  calibrate  both  the  compression  and  shear  response  of 
the  gage,  showed  the  need  for  a  better  understanding  of  several  effects 
not  originally  envisioned.  Some  of  the  difficulties  encountered  in 
interpreting  the  results  are  specific  to  the  shear  gage  usage,  but  many 
of  the  difficulties  are  more  fundamental.  For  example,  the  disagreement 
between  the  &22  va^ue  obtained  from  the  shock  work^  and  the  ultrasonic 
studies^  represents  a  major  uncertainty  in  defining  the  optimal  shear 
direction.  Also,  the  value  of  C^  needs  to  be  established  accurately. 

The  main  findings  of  the  present  work  are  summarized  below: 

(1)  The  experimental  data  on  163°-rotated  Y-cut  LiNbO-j 
confirm  the  large  sensitivity  to  shear  relative  to 
compression.  There  does  not  appear  to  be  a  cross-axis 
sensitivity. 

(2)  The  current  output  from  a  165 . 5°-rotated  Y-cut  gage  is 
considerably  higher  than  the  output  from  a  163°  Y-cut 
gage  for  the  same  compression  loading.  This  result  is 
in  agreement  with  the  theory. 

(3)  Because  of  the  large  shear  sensitivity  of  the  gage, 
proper  interpretation  of  the  experimental  data  requires 
that  the  effects  of  both  tilt  and  anisotropy  be  included 
in  the  analysis.  An  elastic-dielectric  analysis  has 
been  developed  for  incorporating  these  effects  in  a 
quantitative  manner  in  analyzing  the  data;  including 
these  effects  leads  to  significant  changes. 

(4)  Experimental  measurements  of  the  shear  contribution  to 
polarization  agree  well  with  theoretical  predictions  to 
within  3  percent. 

(5)  The  quantitative  agreement  between  theory  and  experiment 
is  not  good  for  compression  loading.  The  inclusion  of 
tilt  effects  makes  a  significant  difference,  but 
significant  discrepancies  (10  to  50  percent)  remain. 
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(6)  The  theoretical  developments  presented  in  this  study 
have  pointed  out  the  complexity  in  analyzing  the 
response  of  rotated-cut  crystals.  In  particular,  it  is 
important  to  develop  a  quantitative  understanding  of  the 
influence  of  the  boundary  conditions  on  the  electrical 
fields,  that  is,  evaluate  the  two-dimensional  electro¬ 
static  effects. 

(7)  Large  electromechanical  coupling  effects  (current 
ramping)  were  observed  in  the  experiments,  but  these 
have  not  been  modeled.  Part  of  the  discrepancy 
indicated  above  could  be  caused  by  electromechanical 
effects.  Future  work  should  incorporate  these  effects 
in  the  analysis. 

(8)  We  have  developed  a  theoretical  analysis  to  incorporate 
electromechanical  coupling  under  combined  compression 
and  shear  loading.  Numerical  procedures  for  analyzing 
experimental  data,  similar  to  the  procedures  in 
Reference  15,  need  to  be  developed. 

(9)  Some  of  the  inconsistencies  in  the  experimental  data 
have  not  been  explained.  Theoretical  calculations  for 
the  163°  Y-cut  orientation  indicate  that  small  errors  in 
crystal  orientations  (±  0.2°)  can  lead  to  large  errors 
in  the  current  output  for  compressive  loading.  This 
result  suggests  the  need  for  an  accurate  knowledge  of 
the  gage  orientation. 

(10)  The  extreme  sensitivity  of  the  gage  output  to  crystal 
orientation  for  compression  loading  suggests  that  the 
gage  orientation  that  is  used  should  be  a  compromise 
between  minimizing  the  compression  output  and  accommo¬ 
dating  the  practical  constraints  on  crystal  orientation 
tolerances.  Fortunately,  the  piezoelectric  constant 
controlling  the  shear  output  is  largely  insensitive  to 
changes  in  crystal  orientation. 

(11)  The  optimal  shear  direction  is  highly  sensitive  to  the 
value  of  the  piezoelectric  constant  e^o*  Analysis  of 
Graham's  impact  data0  gives  an  e22  value  of  2.82  C/mz  in 
contrast  to  the  ultrasonic  value  of  2.43  C/m  .  This 
discrepancy  needs  to  be  resolved.  For  example,  is  it 
reasonable  to  neglect  iE  in  equation  (2-24)  in  analyzing 
the  impact  data?  The  tilt  contribution  in  Graham's 
experiments  is  not  known.  Based  on  our  work,  we 
recommend  using  an  e22  value  of  2.43  C/m  . 

(12)  Because  the  compression  response  of  the  ^63°  Y-cut 
orientation  depends  strongly  upon  the  value,  this 
constant  needs  to  be  established  accurately. 
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These  findings  from  the  present  work  have  resulted  from  a  detailed 
examination  of  the  shear  gage  concept,  which  was  the  main  objective  of 
this  study.  We  believe  that  the  development  of  a  piezoelectric  shear 
stress  gage  is  a  realistic  undertaking.  However,  this  development  will 
require  additional  research  to  understand  the  piezoelectric  response  in 
detail.  Some  specific  recommendations  for  further  work,  including  the 
development  of  a  gage  package,  are  discussed  below. 

We  believe  that  a  basic  necessity  is  to  better  understand  and  model 
the  compression  response  of  the  gage.  This  assertion  may  be  questioned 
because  the  shear  response  has  been  modeled  reasonably  well.  The  main 
reason  for  developing  a  better  understanding  of  the  compression  response 
is  that  in  most  applications  the  gage  will  be  subjected  simultaneously 
to  a  compression  and  a  shear  stress.  Even  if  the  compression  stress  is 
independently  known,  the  determination  of  the  shear  stress  from  the  gage 
output  requires  that  the  sensitivity  to  compression  stress  should  be 
both  small  and  known  accurately.  This  requirement  is  necessitated  by 
two  factors.  First,  in  most  applications  the  compressive  stresses  are 
much  larger  than  shear  stress.  Second,  the  compressive  response  of  the 
gage  is  very  sensitive  to  the  orientation  of  the  crystal. 

The  general  approach  to  studying  the  compression  response  of  the 
gage  should  begin  by  developing  a  quantitative  solution  to  the  two- 
dimensional  electrostatic  problem  presented  in  Section  2.2.  The  highest 
priority  should  be  given  to  this  assessment  of  the  effects  of  the 
boundary  conditions  because  it  is  necessary  for  evaluating  the  complete 
response  of  rotated-cut  crystals  including  the  current  ig  in  Eq.  (2-24). 

This  analysis  should  be  verified  by  modeling  both  shunted  and 
shorted  gage  data  for  identical  impact  conditions  and  comparing  these 
predictions  with  laboratory  measurements  on  such  gages.  Additionally, 
numerical  methods  for  incorporating  into  the  data  analysis  the  electro- 
mechanial  coupling  effects  described  in  Section  2.3  should  be  developed. 

If  these  tasks  can  be  successfully  completed,  then  further 
laboratory  impact  experiments  should  be  conducted  to  determine  the  e22 
constant  very  accurately.  These  experiments  need  to  be  of  higher 
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precision  chan  Che  presenc  experiments  and  should  IncorporaCe  Che 
following  changes:  shunCed  or  shorced  gages  should  be  used  depending  on 
Che  resulCs  of  Che  elecCrosCaCic  analyses;  compression-only  experimenCs 
wichouC  a  buffer  place  should  be  performed;  Che  impaccor  maCerlal  should 
be  more  nearly  elascic;  improved  measuremenCs  of  impacC  tile  should  be 
obcained;  cryscal  orienCaCion  should  be  deCermined  accuraCely  by  several 
independenc  measuremenCs.  After  Che  ^2  constant  is  deCermined,  Che 
recommendaCion  indicaCed  in  iCem  (10)  above  should  be  performed.  The 
crysCal  orienCaCion  Chac  is  selecCed  should  Chen  be  calibraCed  carefully 
under  compressive  loading. 

Regarding  Che  gage  package  developmenc,  an  imporCanc  firsc  step  is 
Co  develop  analycic  or  numerical  mechods  for  approximaCing  Che  scress 
and  scrain  disCribuClons  expecCed  in  Che  gage  in  some  Cypical  applica- 
cions.  Because  of  Che  anisocroplc  response  of  Che  cryscal,  chis  is  an 
involved  undertaking.*  Once  Che  mechanical  quanCiCies  are  known,  Che 
polarization  and  elecCric  field  disCribuClons  can  be  deCermined.  Just 
as  for  studying  the  compression  response  of  LiNbO^  in  the  laboratory,  in 
developing  a  field  gage  package,  it  is  important  Co  develop  a  good 
understanding  of  the  electrostatic  response  of  the  gage  for  the 
appropriate  boundary  conditions  Co  evaluate  deviations  from  the 
idealized  response. 

We  cannot  overemphasize  the  importance  of  identifying  and  under¬ 
standing  all  the  various  contributions  to  the  gage  output  before  using 
the  gage.  Without  a  detailed  understanding  of  these  contributions,  it 
will  be  nearly  impossible  to  analyze  the  gage  data. 

Calculation  of  the  mechanical  quantities  indicated  above  will  also 
be  useful  in  evaluating  the  inclusion  problem:  How  is  the  mechanical 


The  development  of  a  numerical  method  to  address  this  problem  would 
also  be  useful  in  Incorporating  the  crystal  anisotropy  in  the  tilt 
analysis • 


state  of  the  sensor  related  to  the  measurement  of  interest?  The  answer 
to  this  question  is  central  to  all  gage  measurements* 

Finally,  the  development  of  a  piezoelectric  shear  gage  is  dependent 
upon  extending  our  basic  understanding  of  piezoelectric  phenomenon  as 
indicated  in  this  section.  Hence,  it  is  difficult  to  project  accurately 
the  likelihood  of  success  or  time  duration  for  completion  of  an  effort 
to  develop  a  field  gage  package. 
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APPENDIX  A 


SYNOPSIS  OF  THE  FEASIBILITY  STUDY* 

A.l  MOTIVATION  AND  OBJECTIVES 

The  determination  of  dynamic  stresses  and  loads  is  fundamental  to 

much  of  DNA  field  testing.  The  need  for  these  measurements  has  led  to 

the  development  and  use  of  many  different  types  of  stress,  particle 

27  28 

velocity,  acceleration,  and  displacement  gages.  ’  Despite  the  large 
variety  and  quantity  of  existing  dynamic  measurements,  field  measure¬ 
ments  of  shear  stress  (or  loads)  are  lacking.  The  inability  to  make 
dynamic  shear  measurements  is  an  important  shortcoming  because  shear 
measurements  are  needed  for  determining  strength  properties  of  both 
soils  and  structures  in  underground  tests. 

The  need  for  shear  measurements  is  well  recognized  by  most  workers 
involved  in  dynamic  measurements.  However,  this  development  has  been 
lacking  due  to  the  complexity  of  the  problem:  The  desired  shear  stress 
gage  must  be  usable  under  complex  loading  conditions,  and  suitable 
methods  are  needed  to  calibrate  the  gage  to  well-defined  shear  stresses 
Most  of  the  previous  field  gage  techniques  have  been  extended  from 
laboratory  concepts  and  measurements,  but  laboratory  studies  are  also 
lacking  in  dynamic  shear  measurements. 

The  objective  of  our  work  was  thus  to  examine  the  feasibility  of 
developing  a  piezoelectric  shear  stress  gage  for  use  in  DNA  field 
tests.  A  combined  analytic  and  laboratory  experimental  effort  was 


The  Introduction  and  Summary  section  of  our  feasibility  work  described 
in  the  DNA  Report  4870F  (December  L978)  is  presented  here.  To  aid  the 
reader,  this  Appendix  has  its  own  set  of  references.  Because  the 
report  was  written  more  than  five  years  ago,  not  all  the  statements 
reflect  the  author's  current  thoughts. 
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undertaken  to  meet  this  objective.  The  bases  for  this  work  are  recent 
developments  at  SRI  relating  to  the  study  of  dynamic  shear  properties  of 
solids  *  and  the  use  of  piezoelectric  gages  in  studying  dynamic 

1  A  1  »\  A  I 

compressive  stresses.  * 

A. 2  BACKGROUND  AND  APPROACH 

Piezoelectricity  and  the  use  of  piezoelectric  transducers  and 

IQ  OC  3/: 

devices  in  acoustic  applications  is  a  major  field  of  study.  *  * 

However,  the  use  of  piezoelectric  transducers  for  measuring  large  stress 
amplitudes  under  dynamic  loading  is  a  more  recent  and  specialized 
topic.  Studies  during  the  past  decade  have  led  to  the  development  of 
stress  transducers  for  studying  dynamic  compressive  stresses  up  to  tens 
of  kilobars  in  materials  subjected  to  impact,  explosive,  and  radiation 
loads. 10,13,31-34  Studies  by  Graham  and  co-workers  have  clarified 
piezoelectric  response  at  high  stresses  and  established  bounds  on  the 
use  of  stress  gages.  *  '  The  formulation  of  nonlinear  piezoelectric 

constitutive  relations  has  also  received  increased  attention.  *  In 
many  laboratory  and  field  measurements,  a-quartz  is  used  as  the  gage 
material,  but  recent  studies  have  also  been  conducted  using  lithium 
niobate  (LiNhC^).**  For  low  stresses  (below  10  kbar),  the  larger 
electrical  output  of  LiNbO^  is  advantageous. 

Piezoelectric  gages  are  used  in  two  modes.  The  current  or  short- 
circuit  mode  measures  fast  rise  time,  short-duration  stress  pulses,  and 
the  useful  recording  time  of  the  gage  is  the  wave  transit  time  through 
the  gage.  Laboratory  shock  wave  experiments  with  zero  lateral  strains 
commonly  use  the  current  mode.  In  the  charge  or  open-circuit  mode,  the 
gage  is  used  to  record  slower  rise  time,  long-duration  (millisecond) 
pulses.  This  second  mode  is  more  commonly  used  under  field  conditions 
for  recording  stress  pulses  with  wavelengths  much  greater  than  the  gage 
thickness.  That  is,  the  gage  in  the  charge  mode  acts  like  a  static 
transducer  in  equilibrium  with  the  surrounding  material  stresses.  No 
fundamental  differences  exist  between  these  modes  with  regard  to  the 
piezoelectric  response,  and  the  Laboratory  results  are  applicable  to 
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field  usage.  The  differences  in  these  modes  are  operational  and 
reflect  the  mechanical  boundary  conditions  most  suited  for  using  the 
gages  in  different  applications.  Further  discussion  of  these  two 
operational  modes  is  given  in  Reference  31. 

Despite  the  many  laboratory  and  field  developments  in  the  use  of 
the  above-cited  modes,  measurements  to  date  have  been  performed  exclu¬ 
sively  for  compressive  stresses.  In  the  following  paragraphs  we  discuss 
our  approach  for  developing  a  piezoelectric  shear  stress  gage  for  use  in 
dynamic  loading  conditions  in  the  presence  of  complex  stresses. 

The  two  general  requirements  in  the  development  of  a  shear  gage 
are:  (a)  knowledge  of  a  material  phenomenon  (e.g.,  appropriate 

piezoelectric  response)  relating  shear  stress  to  a  measurable  quantity 
and  (b)  the  ability  to  calibrate  the  gage  for  known  shear  stresses. 

The  first  requirement  is  easily  satisfied,  in  principle,  because  of 

'i  C 

the  existence  of  many  shear  transducers  in  the  field  of  ultrasonics. 
Because  the  compressive  stress  gage  is  an  extension  of  ultrasonic 
concepts,  these  concepts  can  be  explored  for  the  development  of  a  shear 
stress  gage.  There  are,  however,  complicating  factors  that  do  not  allow 
a  simple  extension  of  ultrasonic  concepts.  In  ultrasonics,  pure  shear 
waves  are  commonly  created  by  using  the  converse  piezoelectric  effect. 
Thus,  the  measuring  transducer  is  subjected  to  a  pure  shear  motion.  In 
most  dynamic  loading  situations,  a  complex  stress  state  exists — that  is, 
superposed  compression  and  shear  states.  Furthermore,  because  of  the 
large  stresses,  the  material  in  which  the  gage  is  placed  undergoes 
inelastic  deformations,  and  the  relative  magnitudes  of  the  compression 
and  shear  stresses  are  expected  to  vary  over  the  time  range  of  interest. 
For  the  gage  measurements  to  be  useful,  the  electrical  signal  from  the 
gage  should  be  uniquely  related  only  to  the  shear  stress  of  interest. 
Because  of  the  tensorial  nature  of  piezoelectricity,  this  requirement  is 
not  met  by  most  of  the  ultrasonic  shear  transducers. 

"k 

The  electromechanical  coupling  will  be  different  for  the  two  cases. 


To  satisfy  the  above  requirement  in  the  presence  of  complex 
stresses  (e.g.,  combined  compression  and  shear),  we  developed  an 
analytic  approach  that  consists  of  simultaneously  examining  wave 
propagation  and  piezoelectric  relations  for  materials  of  interest.  The 
analytic  approach  chosen  is  general  and  provides  criteria  for  deter¬ 
mining  the  needed  gage  designs.  We  examine  all  possible  crystal 
orientations  to  obtain  the  optimal  directions  for  piezoelectric  response 

and  to  determine  pure  mode  directions  ("specific  directions")  for  stress 
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wave  propagation.  Once  the  optimal  gage  designs  have  been 
analytically  evaluated,  we  can  conduct  experiments  to  verify  and 
calibrate  these  gage  designs. 

The  second  requirement,  calibration  of  the  gages,  is  met  by 
experimentally  studying  the  gage  response  to  one-dimensional  compression 
and  shear  waves  using  the  recently  developed  IMPS*  experimental 
facility.^  In  this  method,  the  specimens  are  subjected  to  varying  but 
controlled  amounts  of  compression  and  shear  stresses.  By  subjecting  the 
gages  to  pure  compression,  we  can  ensure  that  there  is  not  electrical 
output  from  the  compressive  stress.  Increased  amounts  of  shear  stresses 
can  then  be  superposed  to  calibrate  the  response  to  shear.  Because  this 
is  a  feasibility  study,  the  present  scope  of  the  work  is  intended  to 
experimentally  verify  the  theoretical  concepts.  A.  detailed  calibration 
would  be  performed  after  the  feasibility  of  the  shear  gage  has  been 
established. 

A.  3  SUMMARY 

The  objective  of  the  work  reported  here  was  to  examine  the 
feasibility  of  developing  a  piezoelectric  shear  stress  gage  sensitive 
only  to  shear  loading.  Using  a  combined  analytic  and  experimental 
approach,  we  successfully  demonstrated  the  feasibility  of  developing 
such  a  gage. 

*Internal  Measurement  of  P  and  S  Waves. 

^Work  performed  under  DNA  Contract  DNA001-76-C-0384 . 


A  simplified  but  general  analysis  of  the  mechanical  and  piezo¬ 
electric  response  was  performed  to  provide  criteria  for  a  suitable  shear 
gage.  Mechanical  wave  propagation  analysis  in  anisotropic  electic  media 
showed  that,  in  general,  three  waves  are  propagated:  one  quasi¬ 
longitudinal  and  two  quasi-transverse  waves.  Only  for  "specific 
directions"  are  the  propagated  waves  purely  longitudinal  and/or  purely 
transverse.  To  avoid  mechanical  coupling  of  strains,  the  gage  thickness 
direction  should  be  along  a  specific  direction.  The  piezoelectric 
analysis  showed  that  electrical  polarization  for  the  desired  gage  should 
be  one-dimensional  and  along  the  gage  thickness  direction.  Furthermore, 
this  electrical  polarization  should  be  caused  only  by  shear  loads  and 
not  compressive  loads. 

The  adequacy  of  a  particular  crystal  type  for  meeting  the  above 
criteria  can  be  easily  and  efficiently  assessed  by  numerical  calcu¬ 
lations  of  the  specific  directions  and  the  "piezoelectric  matrix"  for 
all  possible  orientations  about  the  three  crystallographic  axes.  When 
these  calculations  were  performed  for  the  different  orientations  of 
a-quartz  and  lithium  niobate  (LiNbO^),  only  the  163°  Y-cut  LiNb03 
appeared  suitable  for  use  as  a  shear  gage.  For  this  orientation,  the 
shear-to-compression  sensitivity  ratio  for  the  polarization  is  enhanced 
by  more  than  two  orders  of  magnitude.  Furthermore,  this  orientation 
deviates  only  2°  from  a  specific  direction  and  therefore  has  minimal 
mechanical  coupling. 

Impact  experiments  were  conducted  under  combined  compression  and 
shear  to  verify  the  use  of  163°  Y-cut  LINbO^  as  a  shear  stress  gage. 

The  results  of  these  experiments  show  that  the  gage,  as  desired,  had 
negligible  sensitivity  to  compressive  stress  and  a  very  large 
sensitivity  to  shear  stress. 

Further  development  of  gages  for  routine  use  in  laboratory  and 
field  measurements  requires  calibration  experiments  that  can  be 
performed  using  the  same  impact  facility.  In  addition,  field  usage 
requires  development  of  proper  packaging  techniques. 


APPENDIX  B 


WAVE  PROPAGATION  ALONG  NONSPECIFIC  DIRECTIONS 

An  important  aspect  of  the  present  work  involved  elastic  wave 

propagation  along  different  crystal  orientations.  In  this  Appendix,  we 

briefly  describe  and  summarize  analytic  developments  pertinent  to  the 
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analysis  of  our  experimental  results.  Studies  by  Borgnis,  Hearman, 

2ft 

and  Johnson  serve  as  a  basis  for  the  developments  presented  here. 

As  indicated  in  Section  2.2,  the  mechanical  wave  propagation  cannot 

be  uncoupled  from  the  external  electric  circuit  for  a  crystal  such  as 

LiNbOj*  However,  the  initial  jump  can  be  analyzed  using  elastic 

constants  at  constant  electric  displacement  (C^j).  Here,  we  present 

some  general  developments,  and  calculations  of  specific  directions;  then 

we  describe  numerical  calculations  for  analyzing  impact  data  along  non- 

D  F 

specific  directions.  The  relationship  between  Cij  and  C“jj  is  also 
derived  in  this  Appendix. 

The  rotation  transformation  of  coordinate  axes  defined  and  used  in 
this  Appendix  is  not  consistent  with  the  other  parts  of  the  report. 


B.l  GOVERNING  EQUATIONS 

2ft 

We  consider  two  coordinate  systems:  x^  corresponds  to  the 
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crystallographic  system  and  x^  corresponds  to  the  wave  propagation 
system.  Wave  propagation  is  described  by 
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where 


0  =  density 

uj  =»  material  displacement  along  the  xj-direction 
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second-order  elastic  constants  in  the  primed  system. 


r* 

ijkA 

The  transformation  matrix  between  the  unprimed  and  primed  system  is 
given  by 


a  *  e'  •  e  (B-2] 

mn  m  n 

+  * 
where  e'  and  e  are  unit  vectors  along  the  x  and  directions, 
m  n  m  n 

Plane  wave  solutions  for  equation  (B— 1 )  for  wave  propagation  along 

a  direction  b'  are  given  by 


u!  =  A'  f(t  -  b'x'/c’) 
i  i  n  n 


(b-3: 


where 
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A^31  displacement  amplitude  along  x^ 
c'  *  wave  velocity. 


Substitution  of  equation  (B-3)  in  (B-l)  gives  the  following  equations: 


Kk  -  hk  o^'2)  K  ' 0 


(b-4; 


where 


XIk 


=  b'b!  C'  „  . 
j  JL  i  jki. 


(B-5! 


and  6  is  the  Kronecker  delta.  For  nontrivial  solutions  of  equation 
(B-4),  we  have  the  condition 


!4k  '  6ik  ^ 


=  0 


(3-6) 


Equation  (B-6)  has  three  possible  eigenvalues,  each  of  which  correspond 
to  an  eigenvector  A'.  Each  of  the  eigenvalues  and  eigenvectors 

ID 

represent  the  wave  velocity  and  displacement  amplitude  associated  with 
one  of  the  three  waves.  These  displacements  are  always  mutually 


orthogonal  but  can  have  an  arbitrary  orientation  with  respect  to  the 

wave  propagation  direction  b' .  Only  for  specific  directions  (or 

isotropic  materials)  are  the  displacement  directions  either  parallel 

(longitudinal  waves)  or  perpendicular  (shear  waves)  to  the  direction  of 
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wave  propagation.  In  general,  there  is  one  quasi-longitudinal  wave 
and  two  quasi- transverse  waves. 

In  our  experiments  we  are  subjecting  the  gage  to  externally  applied 
compressive  and  shear  stresses.  If  the  gage  output  is  to  be  related  to 
a  particular  stress  component  (shear  stress  in  the  present  case),  then 
coupling  of  the  stresses  within  the  gage  must  not  occur.  For  the  linear 
elastic  behavior  considered  here,  this  coupling  can  arise  only  as  a 
result  of  the  anisotropy  of  the  crystal.  Therefore,  the  crystals  used 
for  gage  development  should  be  oriented  along  the  specific  directions 
for  the  compressive  and  shear  waves  of  interest  in  our  experiments. 

In  the  remainder  of  this  Appendix  we  restrict  our  analysis  to  the 
crystallographic  Y-Z  (X2~X^)  plane  of  LiNbO^  crystals. 


B.2  WAVE  PROPAGATION  IN  THE  Y-Z  PLANE  OF  LiNbOj 

Let  the  direction  of  wave  propagation  be  b'  =  (0,  b2>  b^). 
components  of  of  (B-5)  are  given  by 
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The  eigenvalues  Q 


pc'  )  of  equation  (B-6)  are  given  by 
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The  solution  of  equation  (B-9)  gives  the  two  eigenvalues  Q3  and  Q3- 
Q3  has  an  eigenvector  along  (1,  0,  0);  of  the  three  waves  that  can 
propagate  in  any  particular  direction  in  the  Y-Z  plane,  one  is  always 
pure  shear  wave  with  particle  motion  along  (1,  0,  0).  The  eigenvectors 
corresponding  to  Q2  and  can  be  calculated  using  the  usual 
procedures.  We  have  written  a  program  that  calculates  the  eigenvalues 
and  eigenvectors  for  wave  propagation  in  the  Y-Z  plane  for  all 
orientations.  Wave  velocities  needed  in  our  work  were  calculated  using 
this  program.  This  program  is  also  useful  in  determining  the  effects  of 
small  deviations  from  the  specific  directions. 

In  our  work  we  have  used  the  general  procedure  described  by 
Borgnis^  for  calculating  specific  directions.  A  vector  p  that 
represents  a  specific  direction  is  both  parallel  to  t>'  and  is  an 
eigenvector  of  equation  (B-6).  Hence,  p  is  one  of  the  form  (0,  P2,  P3), 
satisfying  the  relation 


P  x  P  =■  p2b^  -  p3 b£  -  0 


and  its  components  are  given  by 


X22  b2  +  X23  b3  =  P2 


(B-10) 


(B-ll) 


X23  b2  +  X33  b3  “  p3 
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Substituting  equation  ( B— 1 1 )  in  (B-10)  and  using  the  values  of  from 

equation  (B-7),  we  can  write  the  equation  for  specific  directions  in  the 


Y-Z  plane 


b  \3  /b  \2 

b^j  (C33  "  2C44  “  C13)  +  3C14  (  bj  j 


b3 

+b^^C13-Cll+2C24^  "C1 


(B-12) 


The  solution  to  this  cubic  equation  gives  the  specif c  directions.  Using 
the  Cjj  values  from  Warner  et  al.,^  we  calculated  the  three  specific 
directions  in  the  Y-Z  plane  for  LiNbO^  as  follows: 


(0,  165.44°,  75.44°) 
(0,  45.99°,  -44.01°) 
(0,  104.08°,  14.08°) 


B. 3  IMPACT  CALCULATIONS  FOR  ARBITRARY  DIRECTIONS 
28 

Johnson  has  described  an  analytic  method  for  calculating  the 
impact  response  of  anisotropic  solids.  He  was  also  the  first  to  suggest 
the  use  of  a  Y-cut  quartz  crystal  for  generating  large-amplitude 
compression  and  shear  waves  in  other  materials.  In  our  present  work  we 
need  to  analyze  the  impact  response  of  various  rotated-cut  crystals. 
Unlike  Johnson,  we  chose  to  use  a  one-dimensional  wave  propagation 
program  to  perform  the  needed  calculations.  This  numerical  method  is 
briefly  described  below. 

Most  one-dimensional  wave  propagation  programs  are  restricted  to 
uniaxial  strain  problems.  The  COPS  code  was  developed  to  include  both 
compression  and  shear  wave  propagation  in  solids.  However, the  code's 
treatment  of  the  governing  equations  is  not  completely  general,  and  the 
material  is  assumed  to  be  isotropic.  In  adapting  this  program  to  our 
present  needs,  we  had  to  make  two  changes:  we  generalized  the  treatment 
of  the  governing  equations  and  wrote  a  subroutine  to  incorporate  an 


anisotropic  constitutive  relation.  In  conformity  with  the  rest  of  our 
work,  we  are  restricting  our  development  to  small  strains.  Changes  to 
incorporate  finite  strain  for  the  elastic  deformations  considered  here 
are  straightforward. 

v  * 

The  direction  of  wave  propagation  is  along  the  x^-axis.  The 
governing  equations  for  the  generalized  one-dimensional  wave  propagation 
can  be  written  as 
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(B-14) 


(B-15) 


(B-16 ) 


Here,  p  is  density,  v^  is  the  particle  velocity  (u^),  and  £  is  the 
specific  internal  energy.  The  energy  equation  is  not  used  in  the 
remainder  of  our  discussion. 

The  artificial  viscosity  relations  are  generalized  by  defining  the 
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viscous  stress  as  follows  (repeated  index  notation  is  not  used): 

duj  ^  5u' 

’ll  '  *  Mi<4xP  isf1  +  "t  “*;>  tep  <B‘17> 

and  are  coefficients  for  quadratic  and  linear  viscosity, 

respectively;  Ax'  is  the  spatial  increment  in  the  wave  propagation 
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solution.  The  viscous  stresses  are  used  in  the  usual  manner. 


We  chose  this  axis  to  conform  to  the  existing  notation  in  the  COPS 
program.  In  using  the  results  for  our  work,  indices  have  to  be 
changed  appropriately. 


The  constitutive  equation  in  the  crystallographic  system,  assuming 
isentropic  deformations  is  written  as 


Ci jki.  ekJt 


(B-18) 


This  is  the  usual  form  of  the  linear  elastic  relation  and  is  derived  by 


expanding  the  internal  energy  function  about  an  initial  state.  The 


elastic  constants  used  in  our  calculations  assume  a  constant  electric 


displacement:  c“j^. 


Equation  (B-18)  can  be  used  in  the  wave  propagation  program  either 
by  transforming  equation  (B-18)  to  the  primed  system  or  by  transforming 
stresses  and  strains  back  and  forth  between  the  two  systems.  We  found 
the  latter  procedure  to  be  convenient.  The  subroutine  ANELAS,  written 


to  Implement  this  procedure  ,  consists  of  the  following  steps: 


(1)  Transform  the  strains  in  the  x^-system  obtained  from  equations 
(B-14)  and  (B-15)  into  the  crystallographic  (x^)  system. 


(2)  Change  the  strains  into  the  matrix  notation. 


(3)  Evaluate  stresses  using  the  matrix  form  of  equation  (B-18). 


(4)  Convert  the  stresses t into  the  tensor  notation  and  transform 

them  back  into  the  x^-system  for  use  in  the  momentum  equations 
(B-13)  for  the  next  time  increment. 


The  conversion  between  the  matrix  and  tensor  notation  can  be 


avoided  by  developing  a  procedure  for  transformation  using  the  matrix 


notation . 


To  minimize  the  modifications  to  the  wave  propagation  program,  we 
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always  chose  x^-axis  as  the  direction  of  wave  propagation.  This 
restriction  coupled  with  the  IRE  convention^  for  rotated-cut  crystals 
leads  to  the  following  rotation  transformation  matrices. 


This  subroutine  is  valid  for  crystals  of  three  of  the  five  classes  of 
the  trigonal  system  and  for  all  crystals  of  higher  symmetry. 
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Rotation  about  x^-axis:  This  notation  defines  the  rotated  Y-cut 
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crystals.  The  x^-  and  xj-axes  are  taken  to  be  coincident.  The 
transformation  matrix  for  going  from  the  crystallographic  to  the  primed 
system  is  given  by 


aij 


0  cos©  sin9 
0  -sinQ  cos0 
10  0 


(B-19) 


where  9  the  angle  of  rotation  about  the  x^-axis.  For  most  of  our 
calculations,  this  rotation  matrix  was  used  with  (  9  *  163°  or 
165.5°).  The  inverse  transformation  is  given  by  the  transpose  of  the 

i  « 

matrix  in  (B-19).  Thus,  x^  »  a^j  Xj  and  x^  *  a^  x^. 

The  results  of  the  numerical  calculation  can  be  adapted  to  the 
convention  used  in  the  main  text  of  the  report  by  incrementing,  as 
follows:  1*  2,  2-»  3,  3  +  1. 

» 

Rotation  about  the  x^-axls:  The  X2~  and  x^-axes  coincide.  The 
transformation  matrix  for  going  from  the  crystallographic  system  to  the 
primed  system  is  given  by 


aij 


sin9  0  cos0 

cos9  0  -sin9 

0  1  0 


(B-20) 


The  inverse  transformation  is  given  by  the  transpose  of  this  matrix. 

I 

Rotation  about  the  x-j-axis:  The  x^-  and  x-j-axes  coincide.  The 
transformation  matrix  analogous  to  the  previous  matrices  is  given  by 


aij 


cos0  sin9  0 

-sin0  cos9  0 

0  0  1 


(B-21) 


The  inverse  transformation  is  given  by  the  transpose  of  this  matrix. 


Several  impact  situations  were  simulated  to  ensure  that  the 

numerical  procedure  was  correct.  One  particular  calculation  involved 
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the  simulation  of  Johnson's  analytic  solution  for  the  case  of  Y-cut 

quartz  being  impacted  by  X-cut  quartz.  In  Figure  B.l  the  free-surface 

velocity-time  profiles  using  Johnson's  notation  are  compared  with  his 

analytic  solution.  As  expected,  the  numerical  solution  showed  the 

rounding  due  to  the  incorporation  of  artificial  viscosity.  The 

amplitudes  are  identical  even  for  the  small  particle  velocity  component, 
» 

V2 •  The  wave  velocities  are  in  excellent  agreement  when  the  midpoint  in 
the  numerical  solution  is  used. 

B.4  RELATION  BETWEEN  AND  cf j 

The  paper  by  Smith  and  Welsh^  lists  only  the  elastic  constants  at 

F  D 

constant  electric  field  C^j.  Because  we  need  C^j  in  our  calculations, 
the  relation  between  the  two  sets  of  constants  is  derived  below.  Taking 
the  electric  field  and  the  strains  as  the  independent  variables,  we 
write 


‘  “ij  "V  eu> 


(B-22 ) 


Differentiating  with  respect  to  strains  and  keeping  the  displacement 
field  constant,  we  have 


(B-23 ) 


Using  the  Maxwell  relations  -  (do. ,/5E  )  =  fdD  /&€..),  we  obtain 

ij  m  m  ij 


(B-24) 
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Figure  8.1.  Free-surface  velocity  in  Y-cut  quartz  impacted  by  X-cut 
quartz.  (The  particle  velocities  have  been  normalized  with 
respect  to  the  impact  velocity,  and  the  propagation  time 
has  been  divided  by  the  sample  thickness.  For  this 
calculation  the  direction  of  wave  propagation  was  chosen 
to  be  along  the  x^-axis  in  conformity  with  Johnson's 
notation.  Note  the  small  value  of  Vj.) 
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Because  of  the  summation  over  m  in  the  last  term  and  because  the 
permittivity  is  a  tensor,  care  must  be  exercised  in  simplifying  equation 
(B-24).  For  principal  cuts  in  trigonal  crystals,  the  permittivity 
tensor  has  only  diagonal  terms,  and  equation  (B-24)  can  be  simplified 


easily.  Using  the  matrix  notation  and  the  definition  of  the  piezo¬ 


electric  constants,  we  can  write  (no  summation  convention) 

e. 


+  la  18  +  2g  2B  3  g  38 


(B-25) 


'11 


'22 


'33 


where  e*  .  is  permittivity  at  constant  strain.  In  writing  equation 
(B-25),  the  following  relation  involving  partial  derivatives  is  useful: 


(2Z3  f$£> 

'■dy'z  ^dz'x  '•bx'y 


-1 


(8-26) 


We  emphasize  that  equation  (B-25)  was  derived  by  assuming  that  the 
electrical  permittivity  tensor  has  a  diagonal  form.  Adapting  equation 
(B-25)  to  rotated  cuts  requires  care. 


D  F 

In  Table  B-l,  we  present  the  Cjj  values  calculated  from  the  Cij 


values  of  Smith  and  Welsh.  To  check  our  calculations,  we  also 


D  F 

calculated  the  C^j  values  using  the  C£j  values  of  Warner  et  al.  Except 


for  C?7,  our  calculated  C|(.  values  match  those  of  Warner  et  al. 


,D 


'12 


'ij 


*After  equation  (B-24)  is  written  for  a  few  terms,  the  general  relation 
can  be  easily  deduced. 


/  -  • 
}\'y 


Table  B-l.  Elastic  constants  at  constant  displacement. 


(1011  N/m2) 


Smith  and  Welsh  (Ref.  5)  Warner  et  al.  (Ref.  4) 


iD 

'll 

2.183 

2.19 

,0 

12 

0.425 

0.371 

,D 

13 

0.764 

0.76 

,D 

14 

-0.148 

-0.147 

.D 

'33 

2.495 

2.52 

,D 

'44 

0.956 

0.95 

APPENDIX  C 


DETERMINATION  OF  PIEZOELECTRIC  CONSTANTS  FOR  LiNb03 

Graham^  performed  an  impressive  series  of  impact  experiments  on 
three  different  orientations  of  LiNbO^  crystals:  Z-cut,  Y-cut,  and 
36°-rotated  Y-cut  samples.  His  experiments  with  the  Z-cut  LiNbO^ 
provide  a  direct  determination  of  the  e33  constant.  Using  the  measured 
e33  value  and  the  e^  value  cited  in  the  literature,  Graham  also 
calculated  the  &22  an<*  ei5  constants  from  his  Y-cut  and  36°-rotated  Y- 
cut  data.  Table  C-l  lists  the  constants  reported  by  Graham,^  Smith  and 
Welsh, ^  and  Warner  et  al.^  Constants  from  other  studies  are  listed  in 
Table  VI  of  Graham's  paper. 

In  most  of  our  work,  we  used  the  results  of  Smith  and  Welsh'*  except 
for  the  e33  constant  for  which  we  used  Graham's  value. ^  This  was  a 
reasonable  procedure  because  the  differences  in  the  other  two  constants 

I 

are  small.  In  Appendix  D  we  use  these  values  to  calculate  the  e^ 
matrix  for  rotated  cuts  of  interest.  Using  all  the  constants  from 
Graham's  work  gives  163.6°  Y-cut  as  the  optimal  direction  for  a  shear 
gage  rather  than  the  163°  Y-cut  orientation  that  we  used  in  our  work. 
This  difference  borders  on  the  accuracy  with  which  these  crystals  can  be 
cut . 


Because  Graham's  data  were  not  appropriate  for  determining  e33,  he 
used  the  e^^  value  given  by  Smith  and  Welsh.  Because  Is  small 

the  other  constants  are  not  very  sensitive  to  Its  value. 
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o 

Table  C-l.  Second-order  piezoelectric  stress  constants  (C/m  ). 


Author 

e22 

e33 

e15 

e31 

Warner  et  al., 

1967  (Ref.  4) 

2.5 

1.3 

3.7 

0.2 

Smith  and  Welsh, 

1971  (Ref.  5) 

2.43 

1.33 

3.76 

0.23 

Graham,  1977 
(Ref.  6) 

2.37 

1.80 

3.83 

0.23* 

When  analyzing  our  impact  data,  we  realized  that  the  analysis  of 
wave  propagation  in  a  nonspecific  direction  requires  a  general  aniso¬ 
tropic  analysis.  Such  an  analysis  for  one-dimensional  wave  propagation 
is  discussed  in  Appendix  B.  In  analyzing  the  Y-cut  and  36°-rotated 
Y-cut  crystals,  Graham ^  used  a  uniaxial  strain  analysis.  Although  his 
analysis  was  approximate,  the  small  deviation  from  uniaxial  strain  was 
assumed  to  be  insignificant.  We  reanalyzed  the  Y-cut  and  the  36°-rotated 
Y-cut  data  of  Graham  using  the  anisotropic  analysis  presented  in  the  last 
Appendix.  Our  procedure,  described  below,  leads  to  different  values  for 
an<^  ei5  Chan  those  obtained  by  Graham. 

For  Y-cut  and  36<’-rotated  Y-cut  crystals  a  normal  impact  results  in 
a  quasi-longitudinal  and  two  quasi-shear  waves,  not  simply  a 
longitudinal  wave  (uniaxial  strain).  For  describing  the  response  of 

f 

these  crystals,  we  consider  the  following  coordinate  system:  X^  is 

» 

along  the  crystallographic  X^-axis,  X2  is  along  the  gage  thickness 

I 

direction,  and  X^  is  chosen  to  form  a  right-handed  coordinate  system. 


See  footnote  on  previous  page. 


For  plate  impact  loading,  all  the  quantities  vary  only  along  the  X2- 
dtrection  and  the  strain  tensor  has  three  independent  components: 

£k £;2>- and  €23 w 

» 

The  polarization  along  the  gage  thickness  direction,  P2,  can  be 
written  using  the  matrix  notation* 


e'  € '  +  e'  e' 
22  2  24  €4 


(2-23 


There  is  no  contribution  from  the  e2f$  term  because  &2fy  is  zero.  In 

his  approximate  analysis  Graham^  Ignored  the  contribution  of  the  second 
term  in  equation  (2-23).  Although  may  be  small,  the  large  value  of 
the  piezoelectric  constants  requires  that  both  terms  be  considered  in 
analyzing  the  data. 

Instead  of  analyzing  all  of  Graham's  experiments,  we  used  the 
following  averaging  procedure. 

(1)  Because  we  are  not  interested  in  constants  beyond  second  order 
we  used  the  six  experiments  for  each  orientation  that  gave  a 
linear  polarization-strain  curve.  The  polarization/strain 
ratio  was  averaged  over  the  six  experiments  to  account  for 
experimental  scatter.  These  values,  in  Graham's  notation,  are 
as  follows: 


Y-cut 


36  -rotated 


-2.38  C/m 


-4.65  C/m 


where  r)  is  a  measure  of  finite  strain. 

(2)  One  experiment  for  each  cut  that  gave  a  P^/n  value  close  to 
each  of  these  values  was  selected  for  simulation  using  the 


The  difficulty  of  relating  the  current  output  to  the  polarization  for 
a  rotated-cut,  discussed  in  Section  2.2,  is  ignored  here  as  It  was  in 
Section  4. 


anisotropic  analysis.  For  the  Y-cut  crystals,  we  chose 
experiment  No.  Q-1056  that  gave  a  value  of  2.39,  and  for  the 
36°-rotated  Y-cut  crystals,  we  chose  experiment  No.  Q-1095  that 
gave  a  value  of  4.64.  Using  the  average  polarization/strain 
values  given  above,  we  calculated  the  polarization  values  for 
each  of  these  two  experiments.  This  procedure  allows  us  to 
average  over  the  various  experiments  and  keep  the  polarization 
calculations  independent  of  Graham's  strain  values. 

22 

Using  our  modified  COPS  code,  we  numerically  calculated  the 

mechanical  variables  for  the  two  experiments  cited  above.  The  LiNbO^ 

elastic  constants  cJj  used  in  these  calculations  were  derived  from  the 
F  J  e  ^ 

Cj\,  values  reported  by  Smith  and  Welstv  (see  Appendix  B).  The  results 

from  these  calculations,  along  with  Graham's  uniaxial  strain 

calculations,  are  shown  in  Table  C-2.  The  deviations  from  uniaxial 

strain  are  small.  However,  as  noted  earlier,  the  contributions  of  the 

term  to  the  polarization  are  significant. 

» 

Using  the  transformations  for  e^j,  we  can  express  equation  (C-l)  in 
terms  of  e^j  values  for  rotations  about  the  X^-axis. 
algebraic  manipulation,  it  can  be  shown  that 

P2  *  e22  (c3  €2  "  c2s  €4}  +  e33  ^  e2  +  ^  e4> 

+  e31  (c2s  e'2  -  cs2e^)  +  e15  (2c2s  e’2  +  c3 

where  c  ■  cos  9 ,  s  ■  sin  9 ,  and  9  is  the  angle  of  rotation  about  the  X-¬ 
axis.  For  Y-cut  (9  ■  0°)  and  Z-cut  (9  ■  90°)  crystals,  equation  (C-l) 
is  considerably  simpler, 

Y-cut  LiNb03 : 

Z-cut  LiNbOj :  P^ 


e22  62  +  ®15  e4 


®33  €2 


(C-2) 


After  some 


(C-l) 

-  cs2  €’) 


*When  the  C^..  values  from  Warner  et  al.  were  used,  the  calculated 
stresses  ana  strains  were  close,  as  indicated  in  Table  C-2. 
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For  other  orientations,  we  have  to  consider  all  the  terms  in  equation 
(C-l). 


Using  the  strain  values  from  Table  C-2,  the  value  of  e^  =  0.23 

2  c  p  e. 

C/m  from  Smith  and  Welsh,  and  ejj  *  1.80  C/m4  from  Graham's  work,  we 
can  write  the  polarization  for  the  two  orientations  as  follows: 

Y-Cut  (Q-1056) :  2.24  e22  -  0.28  e15  -  5.26  C/m2  (c_3) 

36 ’-Rotated  Y-Cut  (Q-1095):  0.99  +  1.34  e22  +  1.78  e15  -  11.47  C/m2 


As  indicated  earlier,  the  polarization  values  on  the  right  side  of  these 
equations  represent  an  averaging  over  six  experiments  for  each  orienta¬ 
tion.  Thus,  there  are  small  differences  between  these  polarization 
values  and  those  reported  in  Graham's  paper  for  these  two  shots. 

Solving  the  equations  for  e^  and  e22  gives 


e^^  ■  3.76  C/m" 
e«_  *  2.82  C/m 


(C-4) 


The  e^  value  is  identical  to  that  of  Smith  and  Welsh^  and  is  less  than 
2  percent  different  from  the  value  reported  by  Graham.**  However,  the 
e22  value  is  considerably  different  from  the  value  calculated  by  Graham. 

Although  we  have  results  for  both  163°-rotated  and  165.5°-rotated 
Y-cut  orientations,  we  made  no  attempt  to  calculate  e^  or  a  complete 
set  of  constants  using  our  experimental  results  in  conjunction  with 
Graham's  data  because  the  two  sets  of  experiments  used  different  types 
of  gages. 


•  *•  «•  .*  +*  V  V  .■*  V  v  .*  V  •  V  ■  *  ■  -  “  »  •  i  •  A  ,  -  «  *  .  h  •  ■  *  *  .*»  A\Vj.Vk  »  A  4*»  A  ,%  > 

■  .S.'.V-  -*.  X>  .V  JV  ■*-'C'l»\V  -*«  .*.  l\  A  -H.  J.  *s.  ^ 


mam# 


The  analysis  of  Graham's  data0  described  here  suggests  that  the 
following  values  be  used  for  the  second-order  piezoelectric  stress 
constants  of  LiNbO^.* 

&22  m  2*82  C/m2 
e^3  ■  1.80  C/m2 
e15  *  3.76  C/m2 
e31  -  0.23  C/m2 

For  developing  a  shear  stress  gage,  this  new  set  of  constants  gives 
160. 5°-rotated  Y-cut,  rather  than  163°-rotated  Y-cut,  as  the  optimal 
crystal  orientation.  This  point  is  discussed  further  in  Appendix  D. 


*Recall  that  the  complexities  in  analyzing  results  for  rotated-cut 
crystals  discussed  in  Section  2.2  have  been  ignored. 
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APPENDIX  D 


TRANSFORMATION  OF  PIEZOELECTRIC  CONSTANTS 

This  appendix  describes  a  simple  and  convenient  method  for 
determining  the  crystallographic  orientations  best  suited  for 
piezoelectric  gage  applications.  Although  considerable  information 
exists  on  piezoelectric  constants  along  particular  directions,  the 
method  presented  here  is  more  general  and  convenient  for  the  applica¬ 
tions  of  interest. 

The  piezoelectric  stress  coefficients  (e..^)  and  the  piezoelectric 

19  ^ 

strain  coefficients  (d^j^)  are  defined  by 


Pi  "  eijk  ejk 


(D-l) 


Pi  -  dimn 


(D-2) 


where  the  piezoelectric  stress  and  strain  coefficients  are  related  by 


eijk  *  dipq  Cpqjk 


(D-3) 


Because  the  piezoelectric  constants  are  third-rank  tensors,  they  are 
transformed  as 


e;  ■  a,  a .  a,  e 
ijk  ip  jq  kr  pqr 


(D-4) 


For  the  rotated  Y-cut  crystals  (rotation  about  the  crystallographic  X- 
axis)  used  in  our  work,  the  rotation  matrix  can  be  written  as* 


This  matrix  is  compatible  with  the  coordinate  system  used  in  the  main 
text.  It  is  obtained  from  equation  (B-19)  by  incrementing  the  indices 
of  the  coordinates  in  the  manner  given  following  equation  (B-19). 


a  =* 


1  0  0  ‘ 
0  cos9  sin9 

,0  -sin9  cos9 


(D-5) 


For  LiNb03, 


L 

combining  (D-4)  and  (D-5)  gives 
*22  '  c3*22  +  *3*33  +  *2s(*31 

.  2  x  ,  3  2S 

-c  se22  +  (c  -  cs  )el5 


J 


+  2  .l5) 


+  cs  (e 


33 


e31> 


(D-6) 


(D-7) 


where  we  have  used  the  matrix  notation;  c  and  s  denote  cos9  and  sin9, 
respectively. 

To  determine  the  optimal  orientation  for  a  shear  gage  we  need  to 
minimize  the  piezoelectric  constant  corresponding  to  the  normal  strain 
along  that  orientation  and  maximize  the  piezoelectric  constant 
corresponding  to  the  shear  strain  of  interest.  By  using  the  numerical 
method  outlined  here,  we  can  quickly  determine  the  optimal  orientation. 
Essentially,  the  method  consists  of  implementing  the  transformations 
indicated  in  equations  (D-3)  and  (D-4). 

In  defining  the  coordinate  rotations  for  our  work,  we  chose  the 

*5 

Institute  of  Radio  Engineers  (IRE)  conventions.  Instead  of  the  usual 
solid  angles.  Me  numerically  solved  equations  (B-l)  and  (B-2)  for 
coordinate  rotations  between  0°  and  180°  about  the  crystallographic  X-, 
Y-,  or  7,-axes.  The  numerical  program  has  the  following  inputs. 

(1)  The  axis  of  rotation,  denoted  as  1,  2,  or  3  for  the 
crystallographic  X-,  Y-,  or  Z-axes,  respectively. 

(2)  The  crystal  type  and  the  piezoelectric  matrix  (numerical 
values)  in  the  crystallographic  system. 

The  above  input  is  used  to  determine  the  rotated  constants  as 
follows: 

(l)  Change  the  piezoelectric  constants  from  the  matrix  to  tensor 
notation. 
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(2)  Construct  the  transformation  matrices  for  the  rotation  angle. 

(3)  Transform  the  piezoelectric  tensor. 

(4)  Convert  the  transformed  tensor  back  to  the  matrix  form. 

(5)  Print  out  the  angle  of  rotation  and  the  transformed 
piezoelectric  matrix. 

(6)  Repeat  the  above  steps  for  all  desired  angles. 


The  differences  in  converting  e  and  d  values  from  matrices  to 
tensors  and  vice-versa  are  included. 

The  piezoelectric  constant  matrix  for  LiNbC^  has  the  following 

form: 


(D-8) 


As  indicated  in  Appendix  C,  we  initially  used  the  piezoelectric 
constants  cited  by  Smith  and  Welsh'’  except  for  the  e^-j  constant,  which 
was  taken  from  Graham's  shock  work.^  Of  all  the  piezoelectric 
constants,  the  measurement  of  e-j^  from  shock  experiments  is  most 
reliable.  Thus,  the  constants  for  the  ’’old  set”  are 


&22  ~  2*43  C/m^ 

e^  =*  1.80  C/ra^ 
e15  =  3*76  C/ra^ 

e^  »  0.23  C/ra^ 


(D-9) 


The  e£2  an^  e^5  values  are  slightly  different  from  Graham's  values  (see 
Table  C-l). 


Our  examination  of  Graham's  work  indicated  that  some  of  his 

assumptions  may  not  be  valid.  Hence,  we  reanalyzed  his  data  as 

described  Ln  Appendix  G.  The  "new  set"  of  constants  as  calculated  by 

the  improved  analysis  have  one  significant  difference:  the  value  of  e?>> 
9 

is  2.82  C/m  .  The  other  constants  were  the  same  as  indicated  above. 
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As  an  example  of  the  numerical  procedure  outlined  in  this  Appendix, 
we  show  the  piezoelectric  stress  matrix  (e^  values)  for  LiNbO-j  rotated 
about  the  crystallographic  X-axis  in  Tables  D-l  and  D-2.  The  values 
shown  are  for  the  two  sets  of  constants  and  for  orientations  between 
162°  Y-cut  and  164°  Y-cut.  A  scan  of  Table  D-l  shows  that,  for  gage 

f  » 

thickness  along  the  X2-axis,  the  constant  e22  minimized  for  163°  Y- 
cut  orientation.  Small  changes  in  orientation  lead  to  significant 

t  t 

changes  in  e22  values.  The  e2^  value  is  relatively  constant. 

The  new  set  of  constants  (only  e22  is  different)  yield  a 

significantly  different  value  of  e22  f°r  the  163°  Y-cut  orientation. 

The  optimal  direction  using  the  new  set  of  constants  is  the  160.5°  Y-cut 

»  -3  2 

orientation  with  an  e22  value  of  3.65  x  10  C/m  .  These  results 
highlight  the  need  for  an  accurate  determination  of  e22* 

Results  shown  in  the  two  tables  are  useful  in  rapidly  determining 
the  effects  of  crystal  orientation  on  the  gage  response.  For  example, 

»  t 

the  624  constant  is  relatively  insensitive  in  contrast  to  the 
constant. 
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OBLIQUE  IMPACT  OF  ELASTIC  PLATES 


In  this  Appendix  we  describe  the  theoretical  analyses  for 
evaluating  the  effects  of  impact  tilt  (or  misalignment)  on  the  results 
of  our  experiment.  This  analysis  is  not  exact  because  the  LiNbO^  is 
treated  as  an  isotropic  solid.  Nevertheless,  it  provides  considerable 
insight  into  the  effects  of  tilt.  Figure  E.l  shows  the  experimental 
configuration  we  analyzed.  The  notation  used  here  is  not  consistent 
with  that  in  the  rest  of  the  report.  However,  the  notation,  as  defined 
here,  is  self-consistent.  The  X^-system  is  defined  such  that  the  x{- 
axis  is  normal  to  Plates  2  and  3.  Plate  1  has  a  velocity  vQ  along  the 
positive  x[-direction.  The  tilt  angle,  9,  is  zero  for  an  ideal 
impact.*  The  situation  depicted  in  Figure  E.l  is  valid  for  both 
compression-only  and  combined  compression  and  shear  experiments;  0  is  a 
measure  of  the  deviation  from  parallelism  at  the  instant  of  impact. 
Hence,  the  initial  inclination  of  the  plates  to  the  direction  of 
projectile  motion  is  immaterial. 

As  indicated  in  Section  3,  not  only  is  the  tilt  angle  important 
in  analyzing  the  results  but  so  is  the  orientation  of  the  line  of  con¬ 
tact  relative  to  the  X^  axes.  Angle  <t> ^  defined  in  Figure  E.l  determines 
the  orientation  of  this  line.  The  X^-system  shown  in  Figure  E.l  is 


The  angle  9  defined  here  denotes  a  different  physical  quantity  than 
that  defined  in  Figure  2  of  Section  3. 

^The  two  angles  <t>  defined  here  and  in  Figure  6  of  Section  3  both 
denote  the  orientation  of  the  line  of  contact,  but  they  are  given  rela¬ 
tive  to  differing  axes  of  reference. 


OiM 


obtained  by  a  counter-clockwise  rotation  through  <J>  about  the  x[-axis 

f  I 

when  looking  from  the  positive  X^-axis  (Xj^  and  X^  coincide).  The  X£- 
axis  is  parallel  to  the  direction  in  which  the  line  of  contact, 
traverses  Plate  2.  The  Xj-systera  is  important  because  most  of  the 
analysis  is  performed  for  this  system. 


In  our  analyses  we  considered  both  "no  slipping"  and  "complete 
slipping"  at  each  of  the  two  interfaces  (between  Plates  1  and  2  and 
between  Plates  2  and  3).  These  represent  the  two  extremes  for 
conditions  at  these  boundaries,  and  the  extension  to  intermediate 
conditions  is  straightforward. 

On  impact  of  Plate  1  with  Plate  2  two  waves  that  propagate  away 
from  the  interface  are  generated  in  each  of  the  plates,  as  shown  in 
Figure  E.2.  We  have  divided  the  analysis  of  the  propagation  of  the  two 
waves  into  Plates  2  and  3  into  five  parts:  Part  I  analyzes  the  impact 
shown  in  Figure  E.2;  Part  II  analyzes  the  interaction  of  the 
compressional  wave  P(l)  with  the  interface  of  Plate  2  and  Plate  3;  Part 
III  analyzes  the  interaction  of  the  shear  wave  S(2)  with  this  interface; 
in  Part  IV,  the  stresses  and  particle  velocities  from  the  first  three 
parts  are  summed  in  the  Xj-systera;  finally,  in  Part  V,  all  the 
quantities  are  transformed  to  the  Xj-systera. 

Because  of  the  numerous  waves  considered  in  the  analysis,  it  is 
helpful  to  establish  the  nomenclature,  then  to  review  the  general 
solution  procedure. 


E.l 


NOMENCLATURE 


Dpi  ■  P-wave  velocity  in  the  ic^  material 
Dgj  »  S-wave  velocity  in  the  i^  material 


um  ■  Normal  particle  velocity  behind  the  mt^1  wave 
vm  ■  Shear  particle  velocity  behind  the  wave 
p^Dp^  a  Compression  mechanical  impedance  for 
the  I*’*1  material 


‘Pi 


-  PiDsi  *  Shear  mechanical  impedance  for  the  ic^  material 
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Angle  with  the  interface  made  by  the  P-wave  in  the  i 
material 


th 


Cf- 


fVi 

3^  *  Angle  with  the  interface  made  by  the  S-wave  in  the  iL 
material 

•  Normal  stress  behind  the  mc^  wave;  m  is  an  odd  number  for 
m 


P-waves 


»  Poisson's  ratio  for  the  ic^  material 
*  Density  of  the  ic^  material 

x_  *  Shear  stress  behind  the  mc^  wave;  m  is  an  odd  number  for 
m 

P-waves 

Note,  X.  and  u  are  unchanged  by  a  shear  wave,  and  x  and  v  are 
unchanged  by  a  compression  wave.  The  subscripts  with  X.  and  x  always 
identify  the  wave  and  do  not  imply  the  indices  commonly  used  with. the 
stress  tensor. 

E.2  SOLUTION  PROCEDURE 

The  oblique  impact  shown  in  Figures  E.l  and  E.2  is  analyzed  using 

the  attached-shock  approximation.  This  approximation  is  valid  for  small 

tilt  angles  because  the  velocity  of  the  point  of  contact  is  much  larger 

£ 

than  the  wave  velocities  in  the  material 


tan©  sin© 


(E-l) 


Two  other  approximations  are  used  to  simplify  the  mathematical  analy¬ 
sis:^  (1)  the  turning  angle  of  the  impact  interface  is  neglected,  as 


$ 

>■>) 


Even  for  the  largest  tilts  in  our  experiments  (9  =  0.1°),  vc  is 
116  mm/^s  in  contrast  to  a  P-wave  velocity  of  3  mm/ps  in  PMMA. 

^For  the  first  part  of  the  five-part  analysis  (given  above)  we  analyzed 
a  more  general  problem  that  avoids  these  assumptions.  Because  of  the 
complexity  of  the  equations,  numerical  results  were  not  computed. 
However,  the  solution  for  the  more  general  problem  had  no  significant 
differences . 
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shown  In  Figure  E.2,  and  (2)  in  choosing  the  S-wave  velocity,  we  neglect 
the  density  changes  behind  the  P-wave. 


On  impact,  compression  and  shear  waves  are  produced  in  both 
Plates  1  and  2.  The  inclination  of  these  wave  fronts  can  be  written 
using  Snell's  law  as 


sina^ 


sin^ 


sina2  sin02 


sine 


v 

o 


(E-2) 


The  particle  motion  behind  each  wave  is  assumed  to  be  that  shown  in 
Figure  E.2.  The  choice  of  the  particle  motion  directions  is  not 
important  because  these  are  determined  from  the  final  solution. 


Throughout  the  discussion  we  use  a  right-handed  coordinate  system 
with  the  X^-axis  coming  out  of  the  plane  of  Figure  E.2.  The  coordinate 
systems  needed  for  analyzing  the  flow  behind  each  wave  front  are 
obtained  by  a  counter-clockwise  rotation  of  the  Xj-system  about  the  X-¬ 
axis.  Tensile  stresses  and  strains  are  taken  to  be  positive.  The 
notation  used  for  shear  stresses  and  strains  is  compatible  with  this 
sign  convention.  The  simplest  procedure  to  ensure  consistency  of  signs 
is  to  determine  the  signs  using  the  flow  equations  and  the  definition  of 
strain. 


The  solutions  in  the  first  three  parts  of  the  analysis  are 
obtained  by  combining  the  linear  elastic  constitutive  relations  with  the 
appropriate  boundary  conditions  at  each  interface.  In  each  part  there 
are  four  unknowns,*  and  four  boundary  conditions  are  required  to  solve 
the  problem.  Two  of  these  boundary  conditions  are  common  to  all 
situations: 


Although  there  are  several  stresses  and  particle  velocities,  the 
problem  always  reduces  to  four  unknowns  because  of  the  additional 
relations  provided  by  the  jump  conditions  and  the  constitutive 
equations. 
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(1)  The  normal  stress  Is  continuous  across  the  planar 
interface. 

(2)  The  normal  particle  velocity  is  continuous  across 
the  interface. 


The  two  remaining  boundary  conditions  depend  on  the  frictional 
conditions  assumed  for  the  interface: 


A.  Complete  Slipping 

(3)  The  tangential  stress  is  zero  along  the  interface 
in  the  first  plate 

(4)  The  tangential  stress  is  zero  along  the  interface 
in  the  second  plate 

B.  No  Slipping 

(3)  Tangetial  stress  is  continuous  across  the 
interface. 

(4)  Tangential  particle  velocity  is  continuous  across 
the  interface. 

C.  Frictional  Sliding 

(3)  Tangential  stress  is  continuous  across  the 
interface. 

(4)  Tangential  stress  is  related  to  the  magnitude  of 
the  normal  stress  through  a  frictional  law. 

Here,  we  have  analyzed  only  cases  (A)  and  (B). 


Note  that  a  considerable  amount  of  tedious  algebra  is  involved  in 
deriving  the  equations  presented  in  this  Appendix.  For  the  sake  of 
brevity,  these  derivations  are  not  shown.  The  following  relations  were 
helpful  in  the  algebraic  manipulations.  For  each  material,  we  can  write 


(E-3) 


using  the  symbols  defined  in  Subsection  E.l. 


E.3  IMPACT  OF  PLATES  I  AND  2  (PART  I) 

In  this  subsection  we  are  concerned  only  with  the  four  Initial 
waves  generated  by  the  Impact  of  Plate  1  with  Plate  2,  shown  In  Figure 
E.2.  We  developed  ti  j  solution  for  this  portion  of  the  impact  problem 
before  we  had  settled  on  a  consistent  nomenclature.  Thus,  in  this  part 
of  the  analysis,  the  Impact  configuration  we  analyzed  Is  slightly 
different  from  that  shown  In  Figure  E.2.  This  part  of  the  analysis 
treats  the  flyer  plate  being  tilted  so  that  the  four  wave  fronts  In 
Figure  E.2  are  In  the  top  half  of  the  figure  (A  reflection  about  the  X^- 
X^  plane  in  Figure  E.2).  Because  the  final  result  does  not  appear  to 
depend  on  how  the  Impact  configuration  is  depicted  in  the  problem,  we 
have  made  no  effort  to  rederive  the  relations  leading  to  the  results 
presented  below. 


E.3.1  Complete  Slipping 


The  stresses  and  particle  velocities  behind  each  of  the  wave 
fronts  and  the  flyer  plate  velocity  are  resolved  parallel  and 
perpendicular  to  the  interface.  For  this  case  it  is  convenient  to 
express  the  equations  in  terms  of  the  stresses.  Using  a  matrix 
representation,  we  can  write  these  equations  as 


Tangential  Stress  Zero  in  Plate  2 
Tangential  Stress  Zero  in  Plate  1 
Normal  Stress  Continuous 
Normal  Velocity  Continuous 


(E-4) 
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here  [A]  is  the  4x4  matrix 


sin2c(2  2v^  -  1 

1 

o 

o 

cos282  2  -  2v2 

sin2a^  2v  ^  -  1 

1 

0 

0 

cosB^  2  -  2v^ 

L 

cos28  2 

sin28  2 

-cos20 

sin28  ^ 

cosa^ 

sing  2 

cosa  ^ 

sin2p  ^ 

V 

Z 

p2 

Zpl 

Zsl  J 

(E-5) 


The  boundary  conditions  that  go  with  a  particular  equation  are  indicated 
in  equation  (E-4).  Note  that  the  indices  for  stress  and  particle 
velocity  identify  the  wave,  but  the  indices  in  the  components  of  matrix 
[A]  identify  the  plate. 


In  the  coordinate  system  associated  with  each  wave  the  stresses 
and  particle  velocities  are  related  as  follows: 


Plate  2: 


Plate  1: 


l 

x2 

X. 


-Vi 

_Zs2V2 

-Z  ,u, 
pi  3 

ZslV4 


(E-6) 


In  these  equations,  the  particle  velocities  are  to  be  taken  as  absolute 
quantities.  Each  of  the  above  stresses  refers  to  the  stress  immediately 
behind  the  particular  wave. 

The  set  of  equations  represented  in  equation  (E-4)  can  be  solved 
to  determine  the  four  stresses.  These,  in  turn,  can  be  used  to 
determine  the  particle  velocities  from  equation  (E-6). 
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E . 3 . 2  No 


to'' 


Slipping 


As  before,  Che  stresses  and  particle  velocities  are  resolved 
along  the  interface.  Because  the  tangential  stresses  are  not  zero  on 
both  sides,  the  equations  are  not  as  easily  uncoupled.  In  the  matrix 
representation  we  have 


[A] 


Tangential  Stress  Continuous 
Normal  Stress  Continuous 
Tangential  Velocity  Continuous 
Normal  Velocity  Continuous 


where  the  matrix  [A]  takes  the  form 


(E-7) 


sin2a 


cos2Q  2 


si^a^ 


-cos20^ 


cos202 

sin202 

-00320^ 

sin20  ^ 

sina2 

cos02 

sina 

cos01 

Zp2 

Zs2 

Zsl 

cosa2 

sing  2 

cosa^ 

sin0  ^ 

Zp2 

Zs2 

Zpl 

Zsl 

Once  the  stresses  are  known,  the  particle  velocities  can  be  calculated 
from  the  jump  conditions  indicated  in  equation  (E-6). 


E.3.3  Results 


To  gain  some  insight  into  the  effects  of  tilt  in  our  problem, 
consider  the  following  example:* 


PMMA:  pQ  ■  1.185  g/cm  ,  Dp  =  3  mm/ps,  Dg  =  1.51  mm/ps 


LiNbO^:  pQ  a  4.64  g/cm  ,  D  =  6.75  mm/ps,  Dg  =  4.61  mm/ps 


We  assume  the  following  representative  values 


(E-9) 


v  ■  0.15  mm/ps,  0  *  2  x  10--3  radian 


(E-10) 


Although  LiNbO^  is  not  an  isotropic  solid,  we  treated  it  as  such  and 
used  the  wave  velocities  pertinent  to  the  crystal  orientation  of 
interest  in  our  work.  The  results  are  as  follows: 


Stress 


amplete 


-4.775  kbar 
-4.778  kbar 
-0.4  kbar 
0.097  kbar 


-4.776  kbar 
-4.778  kbar 
-0.395  kbar 
0.105  kbar 


(E-ll) 


al  -  2.29°,  =  1.15*,  ct2  -  5.16°,  62 


From  these  results  we  draw  the  following  conclusions: 

(1)  The  solution  for  the  slipping  or  no-slipping  cases  are 
nearly  identical.^  A  significant  amount  of  shear  stress 
can  be  generated  even  for  small  tilt  angles. 


This  does  not  match  the  experimental  configuration  treated  In  our 
work. 

+  in 

Fowles  reported  a  similar  result  earlier. 
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(2)  The  normal  stress  values  are  very  close  to  the  values 
predicted  for  9  ■  0°.  To  a  first  approximation, 
compression  stress  behind  the  P-waves  is  independent  of 
9. 


(3)  To  a  first  approximation,  the  magnitude  of  the  shear 
stress  depends  only  on  9  and  is  independent  of  vQ. 


(4)  At  least  for  the  complete  slipping  case,  the  direction 
of  the  shear  motion  can  be  specified  a  priori  without 
solving  the  entire  problem.  The  directions  indicated  in 
Figure  E.2  are  correct. 


E.4  P-WAVE  INCIDENT  ON  THE  INTERFACE  BETWEEN  PLATES  2  AND  3 
(PART  II) 


The  compression  wave  generated  by  the  impact  considered  in  the 
previous  subsection,  denoted  as  P(l)  in  Figure  E.2,  interacts  with  the 
interface  of  Plates  2  and  3.  The  four  waves  shown  in  Figure  E.3  are 
produced  by  this  interaction.  The  orientation  of  the  normal  to  a 
particular  wave  front  is  again  calculated  using  Snell's  law;  we  found 
that  equation  (E-2)  determines  all  the  angles  in  the  problem  including 
those  considered  in  Part  III  of  the  analysis.  As  before,  we  present  the 
equations  in  the  matrix  representation  without  presenting  all  the 
algebraic  manipulation. 


E.4.1  Complete  Slipping 


The  stresses  and  particle  velocities  are  resolved  along  the 


1  nterface. 


-X.1cos2B2 


Normal  Stress  Continuous 


/D  2  \2 

-X.ljr —  I  sin2a.  Tangential  Stress  Zero  in  Plate  2 

\  p2/  (E— 12 ) 

0  Tangential  Stress  Zero  in  Plate  3 


Normal  Velocity  Continuous 


' >  /vA.\.Vv v  v  v  -Ov  v  •  -  ..  -A»V 


Figure  E.3.  P-wave  interaction  with  the  interface.  [The  normals  to  the  wave 
fronts  are  shown  with  the  wave  fronts  themselves  indicated  by 
the  short  line  segments  intersecting  the  normals.  The  arrow  ahead 
of  the  wave  front  indicates  the  direction  of  wave  propagation.  The 
arrow  behind  the  wave  front  indicates  the  direction  of  particle 
motion.  The  incident  wave  is  the  P(  1 )  wave  shown  in  Figure  E.2.] 


where  the  matrix  [A]  is  given  by 


The  stresses  and  particle  velocities  are  related  by 


Plate  2:  X^  =  ~Zp2ul 
X7  =  "Zp2u7 
t8  =  +Zs2v8 
Plate  3:  Xj  =*  -Zp^u^ 
t6  *  +Zs3v6 

As  before,  equation  (E-12)  is  solved  by  determining  [A]-1  and 
multiplying  the  right  side  of  equation  (E-12)  by  it.  The  particle 
velocities  are  calculated  from  equation  (E-14). 


The  equations  for  this  case  were  written  in  terms  of  particle 
velocities 

r”  — i 

Normal  Velocity  Continuous 
Tangential  Velocity  Continuous 
Normal  Stress  Continuous 


[A] 


u7 

-u^cosct2 

v8 

3 

-u^sirui2 

u5 

+Zp2ulcosZ^2 

v6 

Ds2 

+Zs2u1F~  Sin2a2 
u  P2  -1 

where  [A]  is  given  by 


(E-15) 


Tangential  Stress  Continuous 


[A] 


-cosa^ 

sin02 

-cosa^ 

-sinp^ 

m 

3ina2 

cos02 

-sina^ 

COS0  j 

-Zp2cos282 

Zs2sin2f32 

Z^cosZP-j 

Zs3sln2,i3 

°s2 

Zs2  OT  sin2a2 
p2 

Zs2cos2^ 2 

D  3 

Z  -  sin2<x~ 

s3  D  ,  3 

p3 

-Zs3cos203 

(E-16) 


Multiplying  both  sides  of  Equation  (E-15)  by  [A]-*-,  we  obtain  Uy,  Vg, 

U5,  and  Vg  in  terras  of  u^ .  The  stresses  corresponding  to  these  particle 
velocities  can  be  calculated  using  the  jump  conditions  in  equation 
(E-14). 


E.5  S-WAVE  INCIDENT  ON  THE  INTERFACE  BETWEEN  PLATES  2  AND  3 
(PART  III) 

Figure  E.4  shows  the  reflection  and  transmission  of  the  shear 
wave,  denoted  as  S(2)  in  Subsection  E.3,  from  the  interface  between 
Plates  2  and  3.  In  accordance  with  the  approximations  set  forth  in 
Subsection  E.2,  the  four  waves  generated  at  the  interface  by  S(2)  are 
taken  to  have  the  same  propagation  direction  as  the  waves  generated  by 
P(l)  shown  in  Figure  E.3.  This  procedure  simplifies  the  derivation  of 


S(2) 

Incident  Wave 


Figure  E.4. 


S-wave  interaction  with  the  interface.  [The  normals  to  the  wave 
fronts  are  shown  with  the  wave  fronts  themselves  indicated  by  the 
short  line  segments  intersecting  the  normals.  The  arrow  ahead  of 
the  wave  front  indicates  the  direction  of  wave  propagation.  The 
arrow  behind  the  wave  front  indicates  the  direction  of  particle 
motion.  The  incident  wave  is  the  S(2)  wave  shown  in  Figure  E.2.] 


v.v  >v  v  »> 


>  .> 


the  equations.  The  angles  made  by  the  wave  front  normals  will  be  the 
same  as  in  Subsection  E.4. 

Except  for  the  contributions  from  the  incident  wave,  the  other 
terms  can  be  written  directly  from  Subsection  E.4  by  incrementing  the 
wave  identification  numbers  by  4.  Because  the  plates  are  the  same  as 
those  in  subsection  E.4,  the  [A]  matrices  will  also  be  the  same. 

E.5.1  Complete  Slipping 

The  four  equations  for  this  case  are 

Normal  Stress  Continuous 
Tangential  Stress  Zero  in  Plate  2 

(E-17) 

Tangential  Stress  Zero  in  Plate  3 
Normal  Velocity  Continuous 

where  the  A  matrix  is  given  by  equation  (E-13). 

The  jump  conditions  are  given  by  substituting 


[A] 


—  — 

—  — 1 

X11 

T2sin2S2 

T12 

-t2cos282 

X 

0 

9 

X2 

Tio 

-  2—  sinp2 
s2 

Zs2  V2 


(E-18) 


for  the  relation  involving  in  equation  (E-14).  The  wave 
identification  subscripts  for  the  other  four  equations  are  to  be 
incremented  by  4  to  conform  to  Figure  E.4. 


E.5.2  No  Slipping 


The  equations  for  this  case  are 


U11 

-v2sin82 

v12 

a 

v2cos8, 

u9 

Zs2V2sin2e2 

v10 

-Zs2v2cos2®2 

_  _ 

_  M 

Normal  Velocity  Continuous 
Tangential  Velocity  Continous 
Normal  Stress  Continuous 
Tangential  Stress  Continuous 


(E-19) 


where  [A]  is  given  by  equation  (E-16).  The  stresses  are  determined  from 
the  particle  velocities  using  equation  (E-18)  and  the  procedure 
indicated  above. 


E.6  STRESSES  AND  PARTICLE  VELOCITIES  IN  THE  Xj-SYSTEM  (PART  IV) 

In  this  part,  the  stresses  and  particle  velocities  from  the 
previous  parts  are  resolved  and  added  in  the  X^-system.  Because  we  are 
interested  in  the  gage  response,  we  restrict  our  calculations  to  Plate 
3.  At  the  end  of  Part  I  of  this  analysis,  the  stresses  and  particle 
velocities  in  Plate  3  are  zero.  Hence,  in  calculating  the  response  of 
Plate  3  we  are  concerned  only  with  Parts  II  and  III  of  this  analysis. 

As  before,  only  the  results  are  presented. 


We  use  <y  to  denote  stresses  and  wm  to  denote  particle 
velocities  in  the  X^-system. 


E.6.1  Part  II  Contributions 


The  P(5)  wave  is  resolved  in  the  X^-system  using  a  clockwise 
rotation  through  about  the  Xj-axis.  The  rotation  matrix  is 


P(5)  Rotation  Matrix  ■ 


cosa^ 

-sina3 

0 

siting 

cosa^ 

0 

0 

0 

1 

(E-20) 


The  S(6)  wave  is  resolved  using  a  clockwise  rotation  through  P3 
about  the  X^-axis.  The  rotation  matrix  has  the  same  form  as  the  one 
above  with  83  is  substituted  for  03. 

The  components  of  the  stresses  and  particle  velocities  existing 
at  the  end  of  Part  II  of  this  analysis  are  obtained  by  summing  up  the 
contributions  from  the  resolved  P(5)  and  S(6)  waves.  Thus,  in  the 
X^-systera  the  stresses  are 


8 


k 


$ 


all  ”  ^5  cosZPj  -  TgsinzPj 


2  2 

ff22  “  X^sin  <13  +  Q3  cos  33)  +  Xgsin203 


(E-21) 


o33  -  q3  x5 

a  12  *  031  ”  -V2  (Q3  ~  1)X^  sin2a3  +  X£cos203 


where 


The  particle  velocities  are  given  by 


w^  ■  u5C0sct3  +  Vgsin{33 
w2  “  Ujsinotj  -  Vgcos03 


(E-22) 


W3  *  0 


In  these  equations  the  u^,  v^,  X^,  and  x^  refer  to  the  coordinate 


system  associated  with  the  icn  wave. 


E.6.2  Part  III  Contributions 


Because  waves  P(9)  and  S(10)  have  the  same  orientation  as  P(5) 
and  S(6),  respectively,  the  rotation  matrices  resolving  each  of  them  in 
the  Xj-system  are  those  used  above.  The  expressions  for  the  components 


of  crjj  and  wm  existing  at  the  end  of  Part  III  of  this  analysis  are 


obtained  by  incrementing  by  4  the  wave  identification  subscripts, 
appearing  on  the  variables  u^,  v^,  X^,  and  x^,  in  Eqs.  (E-21)  and 
(E-22). 


E.6.3  Total  Contributions  Existing  at  the  End  of  Part  III 


The  total  stresses  and  particle  velocities  in  the  X^-system  are 

ou  -  cos2P3(X5  +  \9)  -  sinjPjCxg  +  x1Q) 

2  2 

«  (sin  a3  +  Q^cos  a3>(X5  +  X 9>  +  sin203(x6  +  x1Q) 


O33  ■  Q3(X3  +  X^) 


(E-23 


a,.  -  -  V2 (Q3  “  1)(X,  +  Xq)sin2a,  +  (x,  +  x.0)cos2g. 


12  U21  ^vx3  ‘'^5  -9'— 3  ■  vv6  ■  .10. 

.  _  *  o~.  »  0 


013  ”  a31 


a23  *  a32 


wL  *  cosa3(u5  +  u?)  +  sin03(v6  +  v^) 


w 2  *  sina3(u5  +  V  "  COsP3(v6  +  V10) 


(E-24 


with  Q  as  defined  above. 


E.7  STRESSES  AND  PARTICLE  VELOCITIES  IN  THE  X^-SYSTEM  (PART  V) 

Finally,  the  stresses  and  particle  velocities  are  transformed  to 
the  xi-systetn  indicated  in  Figure  E.l.  To  go  from  the  X^-system  to 
the  X^-system  involves  a  counter-clockwise  rotation  through  the  angle  <J> 
about  the  X^-axis  (looking  at  origin  from  the  positive  X^-axis).  The 
rotation  matrix  is  given  by 


0 

cos<t> 

-sin$ 


(E-25 


The  stresses  and  particle  velocities  in  the  X^-system  are 


^ij  *  aim  ajn  °mn 
wm  ”  amn  wn 

116 


(E-26 

(E-27 


KfvVV 


where  the  amn  and  wn  are  given  in  equations  (E-23)  and  (E-24), 
respectively. 

We  developed  a  program  called  ELPLAT*  to  implement  the  steps 
outlined  in  this  Appendix;  the  program  closely  matches  the  calculation 
sequence  indicated  here.  The  input  to  this  program  consists  of 
specifying  the  densities  and  elastic  wave  velocities  for  each  of  the 
three  plates,  the  flyer  plate  velocity  vQ  along  the  X^-axis,  the  angles 
0  and  <t> ,  and  the  slipping  condition  at  each  interfaces.  Actually,  only 
the  second  interface  condition  is  important  because  the  solution  in 
Part  I  is  independent  of  the  interface  condition. 

The  analysis  presented  here  is  not  exact  for  our  experiments 
because  Plate  3  (LiNbO^)  is  treated  as  an  Isotropic  solid.  Despite  this 
limitation,  an  assessment  of  the  tilt  effects  on  the  particular  compres¬ 
sion  and  shear  stresses  of  interest  can  be  ascertained.  However,  the 
contribution  to  the  other  stresses  cannot  be  determined.  We  have 
performed  preliminary  work  to  extend  our  analysis  to  anisotropic 
plates.  This  will  be  completed  in  a  future  study. 

In  the  theoretical  analysis  presented  in  this  Appendix,  time  is 
not  included  as  a  variable.  Comparisons  with  time-resolved  measurements 
will  require  a  determination  of  the  wave  transit  times  through  Plate  2. 
When  calculating  these  times,  it  is  important  to  account  for  the  wave 
fronts  not  being  parallel  to  the  interfaces. 


*The  author  is  grateful  to  T.  Radzekewicz  for  his  assistance  in  writing 
the  program.  A  copy  of  this  program  is  available  from  the  author. 
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APPENDIX  F 


DETAILS  OF  MEASURED  WAVE  PROFILES 

Details  of  the  voltage-  or  current-time  profiles  from  the  different 
experiments  are  presented  here.  The  experimental  parameters  correspond¬ 
ing  to  these  profiles  are  listed  in  Tables  1  and  2  of  Section  3. 
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Figure  F.l.  Voltage-time  profile  for  experiment  1  (78-2-46). 
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lies  for  experiment  4  (81-2-23). 
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